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1 .  Introduction 


1.1.  Purpose 

;  This  report11  describes  the  findings  of  the  Algorithm  Analysis  Subtask 
group  working  on  the  0.S.  Army  Intelligence  Center  and  School  (USAICS)  Software 
Analysis  and  Management  System  task  (DSAMS)  regarding  geographic  transformation 
algorithms  used  in  four  of  the  intelligence-gathering  systems  under  0SAICS 
cognizance.  In  this  report  a  set  of  parameters  is  developed  to  characterize 
and  catalogue  Intelligence  system  algorithms  in  four  specific  systems.  Indivi¬ 
dual  algorithms  are  analyzed  to  determine  whether  they  are  performing  their 
functions  properly.  Algorithms  that  perform  the  same  function  in  different 
systems  are  compared  to  determine  which  ones  are  best  according  to  various 
criteria. 


The  algorithms  examined  in  this  report  are  taken  from  the  MAGIIC, 
Guardrail,  Trailblazer,  and  BETA  systems.  They  were  chosen  from  the  approxima¬ 
tely  41  deployed  intelligence  systems  for  which  USAICS  is  Combat  Developer 
because  their  documentation  was  quickly  accessible  and  because  they  represented 
a  range  of  algorithm  applications.  Geographic  transformation  (mapping)  algo¬ 
rithms  were  chosen  for  this  report  since  all  four  systems  contain  position 
location/description  functions  and  many  of  their  algorithms  are  unclassified. 

\  \ 

\ 

\ 

1.2.  Background 

Each  of  the  about  41  intelligence  systems  under  DSAICS  cognizance 
employs  several  types  of  algorithms  to  carry  out  its  gathering  and  processing 
of  intelligence  data.  Two  important  types  of  these  algorithms,  geographic 
transformation  and  correlation,  have  been  chosen  for  analysis  during  this  year. 

The  former  translates  grid  zone  locations,  for  example,  from  latitude-longitude 
to  Universal  Transverse  Mercator  (UTM),  while  the  latter  resolves  many  indivi¬ 
dual  sitings  into  militarily  recognizable  targets  based  chiefly  on  standard 
statistical  procedures.  It  is  important  to  develop  a  set  of  parameters  to 
characterize  these  algorithms  so  that  how  they  should  be  catalogued  can  be 


*Two  additional  reports  will  be  submitted  in  FTS2:  a  correlation  analy¬ 
sis  report  and  a  report  on  possible  algorithm  analysis  methodologies. 


determined.  When  these  activities  are  completed,  it  becomes  possible  to  com¬ 
pare  algorithms  that  perform  the  same  function  in  different  systems. 

To  begin  this  process,  the  JPL  Algorithm  Analysis  Sub-task  group  has 
examined  the  geographic  transformation  algorithms  for  four  of  the  41  systems, 
namely  The  Mobile  Army  Ground  Imagery  Interpretation  Center  (MAGIIC),  Guard¬ 
rail,  Trallblazer,  and  Battlefield  Exploitation  and  Target  Acquisition  (BETA). 
These  four  systems  chosen  for  more  detailed  study  represent  several  intelli¬ 
gence  data  analysis  functions.  MAGIIC  is  a  ground-based  analysis  system  to 
assist  in  interpreting  hard-copy  images  from  different  airborne  surveillance 
systems,  including  a  capability  for  computerized  mensuration  on  imagery;  it  can 
also  receive  and  analyze  data  from  Tactical  Electronic  Intelligence  (TEREC) 
collection  systems  and  provide  emitter  location  estimates.  Guardrail  uses 
airborne  sensor  platforms  to  collect  data  on  Direction  Finding  (DF)  emitters; 
extensive  ground-based  software  is  then  used  to  estimate  the  location  of  the 
units,  such  as  command  posts,  associated  with  these  emitters.  Trallblazer  also 
uses  DF  data  to  estimate  emitter  location.  Its  sensor  platforms  are  essential¬ 
ly  fixed  and  ground-based.  BETA  is  a  Test  Bed  program  for  correlating  data 
received  from  several  types  of  sensor  systems  and  making  target  nominations. 
Both  automatic  correlation  and  aggregation  techniques  and  interactive  graphics 
are  used  in  the  operator's  analysis.  These  systems  would  generally  be 
employed  at  Division  or  Corps  level  or  at  an  Air  Force  Tactical  Air  Control 
Element  (TACE)  or  Allied  Tactical  Air  Force  (ATAF) ;  target  nominations  and 
tactical  situation  reports  would  be  available  to  commanders  and  their  staffs 
from  Brigade  through  Echelons  Above  Corps  (EAC). 

DSAICS  has  cognizance  of  a  large  number  of  algorithms  integral  to 
intelligence-gathering  systems  in  various  stages  of  development  and  deployment. 
The  state  of  "deployment”  of  algorithms  in  the  DSAICS  inventory  ranges  from 
that  of  products  of  research  contracts  not  yet  implemented  in  any  system  to 
those  in  fielded  systems  such  as  Trallblazer  or  Guardrail.  In  the  latter 
systems  the  algorithms  are  documented  in  design  documents  (narrative  English 
and  equations),  and/or  in  machine  readable  design  language,  and  in  code.  Often 
not  all  of  these  forms  of  documentation  are  available  for  any  one  system.  For 
reserrch  algorithms  not  yet  implemented  actual  code,  or  even  detailed  flow 
charts,  may  not  be  available;  thus  analysis  must  rely  solely  on  mathematical 
descriptions. 


-2- 


1 

"Algorithm”  means  any  set  of  rules  for  carrying  out  a  single  concep¬ 
tual  operation  on  a  set  of  data,  such  as  transforming  latitude-longitude  co¬ 
ordinates  to  UTM  or  determining  a  position  from  a  number  of  direction  measure¬ 
ments  taken  at  known  points.**  Algorithms  are  often  hierarchical,  lower-level 
algorithms  often  being  used  to  describe  higher-level  algorithms  and  thereby 
illuminating  their  underlying  logical  structure.  Thus,  results  from  one  algo¬ 
rithm  may  be  data  for  another.  USAICS  is  interested  in  algorithms  performing 
intelligence  data  processing  functions  central  to  their  systems'  mission  and 
those  performing  crucial  support  functions,  such  as  geographic  location,  common 
to  a  number  of  systems.  Data  management  or  mathematical  function  algorithms, 
although  vital  to  the  efficient  functioning  of  the  systems,  are  not  being 
treated  in  these  first  algorithm  analyses. 

t.3«  User  Benefits 

These  analyses  can  benefit  users  in  several  ways.  First,  a  catalog  of 
existing  algorithms  will  help  USAICS  avoid  having  algorithms  redeveloped  for 
new  systems  from  first  principles.  Second,  analysis  of  individual  algorithms 
may,  in  a  few  cases,  identify  deficiencies  worth  correcting  on  the  next  system 
revision.  Third,  and  most  important,  the  comparison  of  algorithms  performing 
the  same  function  in  different  systems  can  lead  to  identifying  guidelines  for 
developing  and/or  selecting  algorithms  to  include  in  new  and  revised  systems. 

Selected  algorithms  from  the  systems  studied  will  begin  to  form  a  library  of 
intelligence  algorithms  with  associated  computer  subroutines  that  will  be 
analogous  to  the  Collected  Algorithms  of  the  Association  for  Computing  Mach¬ 
inery  (ACM).  The  creation  of  such  a  library  is  in  the  spirit  of  Ada+,  the 
Department  of  Defense  language  for  embedded  systems,  and  Ada's  environment. 


**These  conceptual  models  should  be  describable,  although  their  tech¬ 
nical  implementation  is  often  significantly  more  complicated  to  present. 

’''Ada  is  a  trademark  of  the  Department  of  Defense 
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2.  Analyzing  the  Algorithms 


2.1.  Early  Steps 

Since  the  Location  and  Movement  Analysis  System  (LAMAS)  system  docu¬ 
mentation  was  available  first,  our  early  analysis  efforts  were  directed  to  that 
system.  A  preliminary  analysis  of  a  Shortest  Path  Algorithm  was  done  anu 
modeled  in  Pascal  as  an  approach  to  standardizing  representations.  This  algo¬ 
rithm  was  a  variant  of  Dijkstra's  Shortest  Path  Algorithm. 

Later  the  sponsor  decided  that  our  first  emphasis  should  be  on  the 
coordinate  conversion  algorithms  of  the  MAGIIC,  Guardrail,  Trailblazer  and  BETA 
systems.  Three  approaches  to  this  analysis  were  tried  and  evaluated.  The 
MAGIIC  system  has  been  hierarchically  analyzed  for  the  interrelationship  of  the 
algorithms.  The  spheroid  models  of  the  Earth's  oblateness  have  been  examined 
for  all  the  systems.  The  grid  zone  generation  algorithms  have  been  compared 
across  all  four  systems. 

2.2.  Learning  Military  Mapping 

To  analyze  the  first  type  of  algorithm  required  learning  the  military 
grid  system.  This  discussion  identifies  the  various  map  projections  and  mili¬ 
tary  grid  reference  systems  examined  and  how  they  are  interrelated.  The  scope 
of  this  discussion  is  limited  to  only  those  map  projections  and  grid  reference 
systems  pertinent  to  the  MAGIIC,  BETA,  Guardrail,  and  Trailblazer  systems. 

The  map  projections  discussed  are  the  Transverse  Mercator,  Polar 
Stereographic,  Lambert  Conformal  Conic,  and  the  Gnomonic.  The  grid  reference 
systems  used  are  the  Universal  Transverse  Mercator  (UTM),  Universal  Polar 
Stereographic  (UPS),  and  Military  Grid  Reference  System  (MGR).  The  selection 
of  a  map  projection  is  based  on  the  properties  it  preserves  in  the  transforma¬ 
tion  from  a  three-dimensional  spheroid  to  a  two-dimensional  plane.  These 
properties  include  orthogonality  of  latitude  and  longitude,  equal  area  repre¬ 
sentation,  distortion  of  shape,  minimal  change  in  scale  faotor  in  either  east- 
west  or  north-south  directions,  and  representation  of  great  circles  by  straight 
lines.  Since  all  map  projections  are  from  a  spheroid  model  of  the  Earth,  the 
parameters  that  the  spheroid  model  use  are  very  important.  Parlous  spheroid 
models  are  used  for  different  portions  of  the  Earth. 


The  selection  of  a  grid  reference  system  depends  on  the  portion  of  the 
Earth  examined  and  the  resolution  desired.  The  UTM  and  UPS  coordinate  systems 
were  adopted  as  standards  by  the  military  .0  minimize  coordination  problems  due 
to  the  proliferation  of  locally-used  grid  reference  systems.  These  coordinate 
systems  are  most  suitable  for  the  representation  of  large  geographic  areas 
(greater  than  9°  in  latitude  and  longitude).  The  MGR  system  provides  greater 
resolution  when  representing  smaller  geographic  areas  (within  100,000-meter  by 
100,000-meter  squares).  The  MGR  system  can  be  overlaid  on  the  UTM  and  UPS 
coordinate  systems  to  eliminate  ambiguity  due  to  repetitions  of  the  100,000- 
meter  square  identifiers.  The  geographic  reference  system  is  simply  given  as  a 
longitude-latitude  pair.  However,  this  reference  system  of  zones  is  cumbersome 
for  representing  locations  in  good  resolution.  Also,  there  is  an  inconsistency 
in  the  form  of  the  coordinates:  some  applications  use  decimal  degree  notation 
while  others  use  clock-like  representations. 

The  UTM  grid  reference  system  is  valid  for  all  longitudes  over  lati¬ 
tudes  between  84°  North  and  80°  South.  This  area  is  divided  into  rectangles  of 
6°  in  longitude  (zones)  by  8°  in  latitude  (bands),  except  for  the  12°  band  from 
72°  to  84°  latitude.  There  are  60  zones  numbered  from  1  through  60  for  the 
zones  from  -180°  to  +180°  longitude.  There  are  22  bands  lettered  C  through  X 
for  the  bands  from  -80°  to  84°  latitude.  There  are  some  subtle  Irregularities 
in  this  pattern  beyond  56°  latitude  between  0°  and  45°  in  longitude  (see  Figure 
2-1).  The  UTM  grid  reference  system  is  based  on  the  Modified  Transverse  Merca¬ 
tor  projection,  but  can  be  mathematically  transformed  for  use  with  other  types 
of  projections. 

The  UPS  grid  reference  system  is  valid  for  the  North  Polar  (+84° 
longitude  to  the  pole)  and  the  South  Polar  regions  (-80°  longitude  to  the 
pole).  These  regions  have  a  grid  zone  number  of  zero  and  consist  only  of  a 
grid  zone  letter  that  is  longitude-dependent.  The  North  Polar  region  grid  zone 
letters  are  Y  (Western  hemisphere)  and  Z  (Eastern  Hemisphere).  The  South  Polar 
regions  are  A  and  B.  The  UPS  grid  reference  system  is  based  on  the  Polar 
Stereographic  projection  (see  Figure  2-2). 

The  MGR  grid  reference  system,  illustrated  for  the  UPS  system  in 
Figure  2-2  and  for  the  UTM  system  in  Figure  2-3,  provides  finer  resolution  than 
the  UTM  or  UPS  grid  reference  systems.  It  identifies  100,000-meter  by  100,000- 
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meter  squares  by  two  letters,  an  Easting  letter  and  a  Northing  letter.  These 
letters  are  sequenced  so  as  to  provide  at  least  18°  separation  between  simi- 
larily-lettered  squares  (within  a  given  spheroid  model  area  -  otherwise  the 
separation  is  9°).  These  lettering  sequences  are  biased  and  restarted  at  the 
boundaries  of  the  underlying  spheroid  models. 

These  MGR  system  letter  designations  may  be  used  without  reference  to 
the  UTM  or  DPS  designations  when  there  is  no  likelihood  of  ambiguity,  otherwise 
the  DTM  or  UPS  designation  is  included.  Positions  within  these  squares  can  be 
interpolated  in  tens  of  meters  and  are  referred  to  as  Easting  and  Northing 
terms  representing  distances  rather  than  degrees. 

The  Transverse  Mercator  projection  transforms  the  Earth's  spheroid 
onto  a  cylinder  secant  to  the  Earth  and  perpendicular  to  its  axis.  This 
projection  is  used  at  latitudes  within  84°  North  and  80°  South.  Scale 
linearity  is  correct  at  the  two  meridians  cut  by  the  cylinder  (6°  apart)  and 
quite  accurate  in  the  band  formed  by  them.  Because  of  the  vertical  linearity 
this  projection  is  particularily  suitable  to  areas  of  interest  in  the  North- 
South  direction.  This  projection  lends  itself  well  to  being  overlaid  with  a 
rectangular  grid  reference  system  (such  as  the  MGR  system). 

The  Polar  Stereographic  projection  transforms  the  Earth's  spheroid 
onto  a  plane  tangential  to  the  Earth  (at  the  pole,  in  our  applications).  This 
projection  is  used  at  latitudes  beyond  84°  North  and  from  80°  South.  Scale 
linearity  decreases  and  equal  area  exaggeration  increases  as  the  distance  from 
the  pole  increases.  Latitude-longitude  orthogonality  is  preserved  at  the 
meridian  crossings.  All  circles  of  latitude  are  concentric,  centered  at  the 
pole.  Thus,  this  projection  is  useful  for  plotting  radio  waves  and  air  naviga¬ 
tion  with  a  compass. 

The  Lambert  Conformal  Conic  projection  transforms  the  Earth's  spheroid 
onto  a  cone  parallel  to  the  Earth's  axis  and  secant  to  the  earth  at  two  lati¬ 
tudes  referred  to  as  the  standard  latitudes.  The  East-West  scale  linearity  is 
correct  at  the  two  standard  latitudes  and  is  relatively  accurate  in  the  band 
between  these  latitudes.  The  projection  preserves  direction  and  shape  quite 
well  within  and  near  the  standard  latitudes.  Hence,  the  Lambert,  Conformal 
Conic  Projection  is  best  suited  to  East-West  measurements  and  is  useful  for  air 
navigation.  Also,  all  meridians  are  straight  and  intersect  at  the  pole.  This 
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projection  is  oost  applicable  to  the  mid-latitude  region  where  the  cone  is 
secant  to  the  Earth. 

The  Gnomonic  projection  transforms  the  Earth's  spheroid  onto  a  plane 
tangent  to  the  Earth’s  at  the  point  of  interest.  This  projection  is  valid  over 
all  latitudes  and  longitudes.  It  has  the  quality  of  representing  all  great 
circle  arcs  on  the  projection  as  straight  lines.  Since  electromagnetic  waves 
travel  the  shortest  distance  route  (great  circle  arc),  the  Gnomonic  projection 
is  ideally  suited  for  the  presentation  of  direction-finding  lines  of  bearing. 

2.3.  Representing  Algorithms  in  Standard  Form 

To  compare  algorithms  across  systems,  all  algorithms  analyzed  must  be 
translated  into  a  standard  format.  Algorithms  in  the  systems  analyzed  had  been 
coded  in  such  diverse  languages  as  assembly  and  structured  FORTRAN  so  that 
translation  into  a  common  Higher-Order  Language  (HOL)  became  essential  to 
searching  for  common  and  diverse  features.  Publication  ALGOL  was  seen  as  an 
attractive  candidate  because  it  has  been  used  in  the  collected  algorithms  of 
the  ACM  for  algorithm  description.  However,  audiences  outside  the  Applied 
Mathematics,  Numerical  Analysis,  and  Computer  Sciences  communities  are  general¬ 
ly  unfamiliar  with  ALGOL;  and  compilers  for  ALGOL  60  or  ALGOL  68  are  not 
readily  available  in  this  country. 

Pascal,  an  ALGOL-like  language,  has  been  chosen  for  the  primary  repre¬ 
sentation  language  for  the  algorithms  because  it  has  many  of  the  properties  of 
ALGOL  (structure,  strong  typing,  etc.),  it  has  become  familiar  to  a  wide  aud¬ 
ience,  and  high  quality  compilers  are  available  on  many  computers  including  the 
Digital  Equipment  Corp.  VAX  and  many  microcomputers  with  CPM  operating  systems. 
The  last  point  is  important  because  the  VAX  will  be  used  by  both  USAICS  and 
JPL,  and  the  microcomputers  are  similar  to  the  word  processors  and  personal 
computers  at  JPL  and  USAICS.  Among  the  features  of  Pascal  that  contribute  to 
its  clarity  are  the  command  structures,  such  as  "if-then-else"  and  "case",  and 
the  user-defined  data  types.  However,  separate  compilations  of  procedures  to 
support  hierarchical  descriptions  of  algorithms  are  an  implementation-dependent 
extension  rather  than  a  basic  feature  of  the  language.  Because  of  this  and 
other  problems  Pascal  provides  at  best  an  interim  solution  to  the  algorithm 
description  problem. 


Ada  offers  a  long-term  solution.  The  Ada  language  avoids  many  of  the 
shortcomings  of  Pascal  and  has  many  additional  features.  A  stronger  reason  for 
using  Ada  is  that  the  Army  is  likely  to  require  all  new  systems  initiated  after 
1984  to  be  programmed  in  it.  While  no  complete  compiler  for  Ada  is  currently 
available,  there  is  an  interpreter  on  the  project  VAX  computer,  although  it  is 
very  slow.  A  compiler  for  an  incomplete  implementation  available  on  Z80-based 
microcomputers  with  CPM  operating  systems  is  also  available.  Although  this 
compiler  is  not  completely  satisfactory  because  of  the  lack  of  user-defined 
types,  it  is  still  useful  for  some  simple  examples  and  for  comparison  with 
Pascal. 
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3.  Characterizing  and  Cataloging  the  Algorithms 

The  set  of  properties  stored  for  algorithms  in  the  database  should 
characterize  an  algorithm  for  military  application  purposes  without  requiring 
that  the  algorithm  itself  be  retrieved  and  examined.  The  algorithm  property 
selection  is  influenced  by  the  following  ways  the  database  will  be  used.  The 
user  may  wish  a  general  summary  of  what  algorithms  are  in  the  database.  A  more 
likely  use  is  to  look  for  algorithms  that  perform  a  specified  function,  such  as 
position  location.  In  another  dimension,  the  user  may  ask  "What  algorithms  do 
we  have  in  MAGIIC?".  Cataloging  properties  should  be  independent  for  efficien¬ 
cy  of  description.  Two  of  the  properties  chosen,  performance  and  robustness, 
are  not  totally  independent.  Requirements  performance  is  interpreted  here  as, 
"does  the  algorithm  do  what  it  says  it  does?".  Lack  of  robustness  is  inter¬ 
preted  here  as  failure  for  special  value  or  failure  because  of  small  errors  in 
input  data.  Range  checking  of  input  variables  is  an  important  contributor  to 
robustness.  This  is  particularly  important  if  the  algorithm  is  to  be  used  for 
more  than  one  system  where  the  calling  programs  can  not  be  expected  to  protect 
it. 

The  Algorithm  Level  Help  File;  (Figure  3-1)  taken  from  the  Acquisition 
and  Database  Entry  Prompting  Tool  (ADEPT)  User's  Guide,  is  examined  and  inter¬ 
preted  as  a  means  for  describing  the  present  classification  scheme.  This 
information  is  provided  for  each  system  in  which  an  algorithm  is  implemented. 
This  set  of  parameters  is  likely  to  change  as  more  experience  is  gained  with 
its  use. 

Examples  of  PSA  reports  are  shown  In  Appendix  7.1. 


Fig.  3-1:  Algorithm  Level  Help  File 


NAME  Is  a  name  description  of  the  algorithms  function. 

SYNONYM  is  the  algorithms  abbreviation  (the  same  as  its  VAX  file 
element  name). 

SOURCE: AUTHOR  is  the  document  from  which  the  algorithm  was  taken  and 
author,  if  known. 

PROCESSING  is  what  JPL  has  done  with  the  algorithm,  e.g.  Pascal  pro¬ 
gram  tested. 

MATH: FIELD  is  what  mathematical  field  the  algorithm  is  based  on,  e.g. 
least  squares. 

ROBUSTNESS  is  a  measure  of  sensitivity  to  values  of  variables  input  to 
the  algorithm,  e.g.  a  transformation  algorithm  may  produce 
an  incorrect  result  with  inputs  of  +  180°  longitude. 

TREE: LEVEL  is  a  rough  Indication  of  location  of  the  algorithm  in  the 
hierarchy  of  algorithms  in  a  system. 

REQUIREMENTS/PERFORMANCE  Since  the  requirements  documents  for  parti¬ 
cular  algorithms  are  generally  not  available,  requirements 
must  be  derived  from  design  documents  or  comments  in  code. 
This  item  describes  bow  well  the  algorithm  meets  these 
requirements . 

REFERENCE  is  a  pointer  to  the  VAX  file  containing  the  representation 
of  the  algorithm  in  standard  form:  Pascal  and  in  some  cases 
Ada.  (These  are  given  in  Appendix  7.3). 
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4.  Algorithms  in  MAGIIC 


The  MAGIIC  System  Lambert  Constant  Generation  algorithm  was  analyzed 
and  modeled  in  Pascal.  This  algorithm  is  required  for  analyzing  and  modeling 
the  Geographic  to  Lambert/Polar  Grid  and  Lambert/Polar  Grid  to  Geographic 
conversions. 

The  MAGIIC  system  coordinate  conversions  were  studied,  and  two  inter¬ 
related  algorithms  were  analyzed  and  modeled  in  Pascal.  These  were  the  Polar 
Grid  to  OPS  and  Northing  and  Easting  to  OTM  conversions.  Several  inconsisten¬ 
cies  have  been  noted  (perhaps  only  because  of  the  technical  writing).  Modeling 
these  interrelated  algorithms  in  Pascal  led  to  the  decision  to  use  a  non¬ 
standard  Pascal  feature  -  the  VAX  Pascal  external  procedure  capability.  This 
was  considered  necessary  to  best  maintain  structural  integrity,  accuracy  and 
clarity  in  the  algorithm  representation. 

4.1*  MAGIIC  Lambert  Constant  Generation  Algorithm 

The  MAGIIC  Lambert  Constant  Generation  Algorithm,  described  in  docu¬ 
ment  CG108100A,  dated  23  October  1978,  paragraph  3.2.119,  page  210,  has  insuf¬ 
ficient  input  parameters  and  a  lack  of  detail  on  setting  the  hemisphere  flags 
to  implement  the  algorithm  in  Pascal  without  making  several  assumptions.  If  we 
assume  the  availability  of  the  underlying  spheroid  parameters,  the  Lambert 
Constant  Generation  Algorithm  performs  satisfactorily.  The  hemisphere  flag 
issue  remains  unsettled:  are  they  north-south  hemispheres,  east-west  hemi¬ 
spheres,  or  both?  -  all  three  selections  make  sense  in  different  contexts. 

This  algorithm  has  been  modeled  in  Pascal  on  the  VAX  computer  for 
uniformity  of  presentation  and  for  comparison  and  analysis.  All  assumptions 
are  included  in  comments  in  this  Pascal  representation  (see  Appendix  7.3). 

4.2.  An  Interrelated  Set  of  MAGIIC  Coordinate  Conversion  Algorithms 

There  are  a  set  of  four  interrelated  coordinate  conversion  algorithms 
for  MAGIIC  that  warrant  a  consolidated  discussion  because  although  they  are 
pair-wise  reciprocal  and  criss-cross  "call"  each  other,  their  input  and  output 
parameters  are  specified  inconsistently.  These  algorithms  viewed  as  reciprocal 
pairs  are: 
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1.  Polar  Grid  to  OPS 

2.  OPS  to  Polar  Grid 

3.  Northing  and  Easting  to  OTM 

4.  OTM  to  Northing  and  Easting 


(paragraph  3.2.118)  reciprocal 
(paragraph  3.2.117)  reciprocal 
(paragraph  3.2.86)  reciprocal 
(paragraph  3.2.85)  reciprocal 


They  can  also  be  viewed  as  "criss-cross"  calling  algorithms  as 

follows : 

1.  Polar  Grid  to  OPS  (paragraph  3.2.118)  each  calls  the 

other 

2.  Northing  and  Easting  to  OTM  (paragraph  3.2.86)  each  calls  the 

other 

3.  OPS  to  Polar  Grid  (paragraph  3.2.117)  each  calls  the 

other 

4.  OTM  to  Northing  and  Easting  (paragraph  3.2.85)  each  calls  the 

other 


Analysis  of  these  algorithms  and  their  interrelationships  has  revealed 
the  following  inconsistencies: 

1.  the  Polar  Grid  to  OPS  algorithm  should  have  an  output  list  con¬ 
sistent  with  the  OPS  to  Polar  Grid  algorithms  input  list  since 
they  are  reciprocal  functions.  The  lists  are  not  consistent; 

2.  the  Northing  and  Easting  to  OTM  algorithm  and  the  OTM  to  Northing 
and  Easting  algorithm,  similarily,  have  inconsistent  output  and 
input  lists; 

3.  in  both  of  the  above  cases  the  same  problems  are  true  when  the 
algorithms  are  used  in  the  "criss-cross  call"  sense. 

These  points  will  be  discussed  in  greater  detail  below  in  the  indepen¬ 
dent  analysis  of  the  four  algorithms. 
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4.3*  MAGIIC  Polar  Grid  to  Universal  Polar  Sterographlc  Algorithm 

The  MAGIIC  Polar  Grid  to  Universal  Polar  Sterographlc  (UPS)  Algorithm, 
described  in  document  CG108100A,  dated  23  October  1978,  paragraph  3.2.118,  page 
209 »  fails  to  produce  the  correct  output  data. 

This  algorithm  first  determines  if  the  UPS  or  UTM  grid  reference 
system  is  applicable  based  on  the  input  grid  zone  designation.  The  algorithm 
that  handles  the  UTM  conversion  is  discussed  as  part  of  the  Northing  and 
Easting  to  Universal  Transverse  Mercator  conversion  algorithm  (paragraph  3.2.86 
of  the  previously  referenced  document). 

In  an  attempt  to  render  this  algorithm  amenable  to  analysis  and  model¬ 
ing  in  Pascal  several  assumptions  have  been  made  (based  on  our  interpretation 
of  the  composite  of  various  inferences  from  paragraph  3.2.86,  .87,  .117,  .118). 
These  assumptions  are: 

1 .  The  Input/output  lists  consist  of  decimal  and  integer  numbers,  and 
alphabetic  letters.  There  are  no  wholesale  conversions  of  the 
input/output  lists  from  and  to  ASCII  representation, 

2.  The  input/output  lists  are  exchanged  with  the  Northing  and  Easting 
to  UTM  algorithm  for  processing,  when  the  input  grid  zone  letter 
is  from  C  thrugh  X  (not  in  either  polar  region).  Otherwise,  this 
algorithm  performs  the  processing  (for  the  polar  regions). 

The  following  discussion  only  treats  the  UPS  coverage  area  where  a 
conversion  from  PS  to  UPS  should  be  made.  This  conversion  is  handled  within 
the  Polar  Grid  to  Universal  Polar  Sterographlc  conversion  algorithm  (paragraph 
3*2.1 18).  This  algorithm  is  defective  in  producing  disallowed  100,000  meter 
squares  letter  pairs. 

The  North  Polar  Area  (i.  84°  N)  and  South  Polar  Area  (>  80°  S)  are 
referenced  using  the  UPS  grid  zone  reference  system  with  the  grid  zone  number 
set  to  zero  followed  by  one  of  the  grid  zone  letters  A,  B,  Y,  or  Z.  Letter 
pairs  represent  the  100,000-meter  squares  which  overlay  the  grid  zone  designa¬ 
tions.  It  is  in  the  lettering  of  these  squares  that  this  algorithm  falls. 


The  descriptions  of  how  to  obtain  the  Easting  letter  Lg  from  the  index 
Ig  and  the  Northing  letter  Ljg  fron  the  index  Ijj  are  merely  statements  that  the 
functions  should  be  performed  without  any  details.  Perhaps  an  elaborated 
description  of  the  details  of  these  functions  would  lead  to  satisfactory  con¬ 
versions,  if  these  details  were  furnished. 

4.4.  MAGIIC  Northing  and  Easting  to  Universal  Transverse  Mercator  Algorithm 

The  MAGIIC  Northing  and  Easting  to  Universal  Transverse  Mercator  (UTM) 
Alorithm  as  described  in  document  CG18100A,  dated  23  October  1978,  paragraph 
3.2.86,  page  156,  appears  to  produce  the  correct  output  data  on  the  UTM  map 
segment  available  for  reference.  It  will  certainly  fail  in  grid  zones  31V, 
3 2X,  34X,  and  36X  and  should  be  further  evaluated. 

This  algorithm  first  determines,  whether  the  UTM  or  UPS  grid  reference 
system  is  applicable  based  on  the  iput  grid  zone  designation.  The  algorithm 
that  handles  the  UPS  conversion  is  discussed  as  part  of  the  Polar  Grid  to 
Universal  Polar  Sterographic  conversion  algorithm  (paragraph  3.2.118  of  the 
previously  referenced  document). 

In  an  attempt  to  render  the  algorithm  amenable  to  analysis  and  to 
model  in  Pascal  several  assumptions  have  been  made  (based  on  our  interpretation 
of  the  composite  of  various  inferences  from  paragraph  3.2.86,  .87,  .117,  .118). 
These  assumptions  are: 

1.  The  input  lists  consist  of  decimal  and  integer  numbers  and  alpha¬ 
betic  letters.  There  is  no  wholesale  conversion  of  the  in  ut 
list  to  ASCII  representation.  The  output  lists  are  entirely  in 
ASCII  representation, 

2.  The  input/output  lists  are  exchanged  with  the  Polar  Grid  to  UPS 
algorithm  for  processing,  when  the  input  grid  zone  letter  is  not 
from  C  through  X  (for  the  polar  regions).  Otherwise,  this  algo¬ 
rithm  performs  the  processing  (for  the  non-polar  regions). 

Only  the  conversion  to  the  UTM  coverage  area,  that  is,  the  actual 
conversion  handled  within  this  Northing  and  Easting  to  Universal  Transverse 
Mercator  conversion  algorithm  (paragraph  3.2.86),  is  discussed  here.  This 
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algorithm  is  defective  in  producing  grid  zone  designation  31V,  which  should  be 
truncated  at  3°  E,  and  the  three  non-existent  grid  zone  designations  32X,  34X, 
and  36X. 


5.  Comparing  the  Algorithms  Across  Systems 

5.1.  Spheroid  Models 

Because  the  Earth  is  not  a  perfect  sphere,  it  is  modeled  as  a  sphe¬ 
roid.  Since  this  underlying  spheroid  model  of  the  Earth's  oblatness  spans  most 
of  the  coordinate  conversion  algorithms,  the  use  of  various  spheroid  models  was 
investigated  for  all  four  systems.  It  was  found  that  twelve  different  spheroid 
models  were  used  among  or  in  the  four  systems  raising  the  possibility  of 
inconsistencies  that  may  hamper  inter-system  communication.  Of  these  spheroid 
models  five  were  used  in  common  by  all  four  systems.  The  remaining  seven 
spheroid  models  were  not  available  in  all  systems,  as  illustrated  in  Table  1. 
Some  of  these  partially-shared  spheroid  models  may  be  equivalent,  but  this 
cannot  be  determined  until  the  code  is  available  for  all  four  systems. 

5.2.  Grid  Zone  Generation 

The  Grid  Zone  Generation  algorithms  were  analyzed  for  the  MAGIIC, 
Guardrail,  Trailblazer,  and  BETA  systems.  The  Guardrail  Grid  Zone  algorithm 
text  description  is  consistent  with  the  Trailblazer  text  description  and  com¬ 
puter  code.  The  MAGIIC  text  description  differs  from  the  Guardrail  and  Trail- 
blazer  descriptions  and  would  tend  to  produce  less  efficient  runtime  code,  but 
would  economize  memory.  All  three  of  these  systems'  algorithms  would  produce 
the  same  flawed  results,  except  BETA  which  fails  badly  at  the  upper  latitudes. 
These  three  versions  of  the  Grid  Zone  Generation  algorithm  have  been  modeled  in 
Pascal. 


The  BETA  Grid  Zone  Generation  algorithm  handles  details  of  grid  zone 
generation  more  completely.  The  grid  zone  number  calculations  handle  input 
longitude  beyond  -180°  and  at  180°  and  beyond  whereas  the  three  other  systems 
would  fail.  The  grid  zone  letter  calculations  handle  the  special  conditions  of 
grid  zone  truncation  for  grid  description  31V  and  the  non-existence  of  grid 
designations  32X,  34X,  and  36X. 

Pascal  Implementations  of  these  algorithms  are  given  in  Appendix  7.3. 
The  underlying  assumptions  are  given  in  the  comments  included  in  the  code. 
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5.3.  MAGIIC  Grid  Zone  Generation  Algorithm 

The  MAGIIC  Grid  Zone  Generation  Algorithm  as  described  in  document 
CG108100A,  dated  23  October  1978,  paragraph  3.2.90,  page  167  has  been  analyzed 
and  found  to  handle  the  following  five  areas  incorrectly: 

1.  the  upper  limit  of  longitude  (180°  E), 

2.  the  latitudes  2  80°  N  (considerably  below  the  upper  limit,  84°  N), 
where  it  provides  erroneous  data, 

3.  the  truncated  grid  zone  31V, 

4.  the  non-existent  grid  zones  32X,  34X,  and  36X, 

5.  the  regions  beyond  the  stated  longitude  and  latitude  limits,  where 
it  fails  catastrophically. 

In  general,  perhaps  due  to  technical  writing,  there  are  many  errors  of 
omission  and/or  commission  where  the  criteria  for  certain  algorithm  parts  are 
left  unstated;  for  example,  a  text  states,  "If  the  latitude  is  84°  north  or 
greater,  or  80°  south  or  greater,  the  grid  zone  number  (sic)  shall  be  3et  to  Y 
or  Z  or  to  A  or  B,  respectively.  The  grid  zone  number  shall  be  set  to  zero." 

5.4.  Guardrail  Grid  Zone  Generation  Algorithm 

The  Guardrail  Grid  Zone  Generation  Algorithm  as  described  in  document 
ESL-TM928,  dated  15  Septemer  1979,  paragraph  16.6.2.1,  page  16-192  fails  in 
five  areas: 

1.  at  longitudes  equal  to  and  beyond  180°  E  and  beyond  180°  W. 

2.  at  latitudes  equal  to  and  beyond  80°  N  and  beyond  80°  S, 

3.  at  the  truncated  grid  zone  3IV, 

4.  at  the  non-existent  grid  zones  32X,  34X,  and  36X, 
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5.  beyond  the  stated  longitude  and  latitude  limits,  where  it  fails 
catastrophicaly. 

5.5.  Trailblazer  Grid  Zone  Generation  Algorithm 

The  Trailblazer  Grid  Zone  Generation  Algorithm  is  described  in  the 
well-commented  ROLM  assembly  language  listings.  This  algorithm,  extracted  from 
the  code  for  the  GP2UH  subprogram  dated  20  February  1981,  fails  in  five  areas: 

1.  at  longitudes  equal  to  and  beyond  180°  E  and  beyond  180°  N, 

2.  at  latitudes  equal  to  and  beyond  8°0  N  and  beyond  80°  A, 

3.  at  the  truncated  grid  zone  3IV, 

4.  at  the  non-existent  grid  zones  32X,  34X,  and  36X, 

5.  beyond  the  stated  longitude  and  latitude  limits. 

Because  the  hierarchical  program  structure  is  not  yet  available  for  Trail- 
blazer,  it  is  possible  that  some  or  all  of  these  problems  are  handled  adequate¬ 
ly  in  higher  levels  of  the  program  structure. 

5.6.  BETA  Grid  Zone  Generation  Algorithm 

The  BETA  Grid  Zone  Generation  Algorithm,  described  in  document  SS22- 
43,  Appendix  IV,  page  11-474  for  the  ADSONU  subprogram,  and  page  II-45P  for  the 
ADSCCM  subprogram,  was  in  Structured  FORTRAN  with  in-line  coding.  The 
"INCLUDE"  subprogram  ZDBPRO  was  missing  so  some  "reasonable"  assumptions  were 
made  about  it. 

This  algorithm  performs  as  specified  and  effectively  handles  the 
following: 

1.  Grid  zone  wrap-around  (longitudes  >180°  W  or  2  180°  E), 


2 


North  and  South  Polar  Regions  (latitudes  2  84°  N  or  >  80°  S), 


3.  Truncating  grid  zone  31V, 

4.  The  non-existent  grid  zones  32X,  34X,  and  36X. 


-22 


Table  5-1:  Inconsistent  Spheroid  Osage 

(X  =  spheroid  model  used  in  the  system) 


Snheroid  Model 

System 

Technics 

BETA 

MAGIIC 

GR 

TB 

TM241 

Clark  1866 

X 

X 

X 

X 

X 

International 

X 

X 

X 

X 

X 

Clark  1880 

X 

X 

X 

X 

X 

Everes  t 

X 

X 

X 

X 

X 

Bessel 

X 

X 

X 

X 

X 

Australian 

X 

X 

X 

X 

Walbeck 

X 

Fisher 

X 

X 

Krasovsky 

X 

X 

World  Geodetic 

X 

X 

Airy 

X 

Malayan 

X 

Reference 

1 

2 

3 

4 

1.  DD2642,  dtd  20  Feb  81,  pg  263 

2.  CG1808100A,  dtd  23  Oct  78,  pg  168 

3.  ESL-TM929,  dtd  15  Sep  79,  Pg  15-158 

4.  TM 32-581 1-022-10-0 
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6.  Discussion  and  Conclusions 


6.1.  Documentation  of  Algorithms 

Only  a  design  document  has  been  available  for  the  MAGIIC  system;  and 
the  code  has  not  been  available,  as  shown  in  Figure  6-1.  Therefore,  some  of 
the  apparent  deficiencies  in  the  algorithms  may  be  due  to  poor  technical  writ¬ 
ing  and  may  not  exist  in  the  code  itself.  For  the  Guardrail  system,  tapes 
containing  code  have  been  available,  but  the  printouts  obtained  to  date  have 
been  largely  unreadable  so  that  suppositions  made  from  the  documentation  could 
not  be  confirmed  from  the  code.  In  the  case  of  Trailblazer,  code  is  available, 
but  the  structured  overview  expected  from  documentation  is  not  available. 

6.2.  Similarity  of  Functions  Across  Systems 

The  functions  performed  by  the  geographic  transformation  algorithms 
are  found  to  be  basically  the  same  across  the  four  systems  examined,  although 
the  functions  are  implemented  in  slightly  different  ways. 

6.3.  Incompleteness  of  Algorithms  for  Global  Applications 

The  MAGIIC,  Guardrail,  and  Trailblazer  transformation  algorithms  do 
not  account  for  all  the  vagaries  of  the  military  grid  system.  Only  the  BETA 
algorithms  account  for  all  regions  and  boundaries.  The  former  systems  may 
ensure  that  "bad"  arguments  are  never  passed  to  these  algorithms,  so  that  no 
anomalies  would  occur  in  overall  system  performance.  However,  to  develop  a 
library  of  algorithms  shared  by  many  systems  requires  that  algorithms  internal¬ 
ly  protect  themselves  from  "bad"  input  data. 

6.4.  Robustness  of  Algorithms 

All  systems  but  BETA  fail  to  check  for  limits  of  latitude  and  longi¬ 
tude.  This  may  arise  from  a  common  tendency  to  focus  attention  on  certain 
areas  of  the  world,  e.g.  Western  Europe.  This  tendency  is  especially  inappro¬ 
priate  when  developing  software  that  may  well  outlive  any  given  political  or 
geographical  constraints. 
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6.5.  Selection/Consolidation  of  Algorithms 


The  BETA  grid  zone  generation  algorithm  is  superior  to  those  in  the 
other  three  systems.  Selection  of  a  spheroid  model  for  the  library  is  not 
possible  on  the  basis  of  the  presently  available  data  and  may  eventually  re¬ 
quire  developing  algorithms  based  on  our  experience  with  many  different 
systems. 
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Fig.  6-1:  Documentation  of  Algorithms 
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7.1.  Database  Entries  and  Products  (PSL/PSA) 

The  three  attached  PSA  Reports  were  found  to  be  useful  to  our  analy¬ 
sis.  The  first  report  is  essentially  the  PSA  Report  representing  our  PSL  input 
data  descriptions.  The  second  is  the  PSA  Data  Activity  Interaction  Matrix.  It 
shows  the  interrelationships  between  the  algorithms  and  their  input  (R)  and 
output  (D)  data  items.  The  third  is  part  of  the  PSA  Structure  Chart  and  is  in 
indented  hierarchy  chart  for  the  algorithms  in  our  PSL  database. 
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This  algorithm  determines  the  appropriate  spheroid  nut fiber 
co  r  respond  iiia  to  the  input  1  a  t  i  tude/ lunSi  lode  be  scanning  a 
specially  constructed  map  database  which  contains  the  relationship 
between  longitude  bands/latitude  strips  and  spheroid  numbers. 


algorithm 

3.2.91  CGI 08100A / P  a  r  t  _ I 


U 

JWG 


abbreviation 
tspe_of_source 
d  a  t  e  _  a  c  a  u  i  [■  e  d 
p  r  o  c  e  s  s  i  n  d  _  d  o  n  e 
m  a  t  h  e  m  a  t  i  c  a  1  _  F  i  e  1  d 
tree-level 

r  e  o  u  i  r  e  m  e  n  t  s  ..  p  e  r  f  o  t-  m  a  n  c  e 


VALUE  : 

Spheroid-Generation 
d  o  c  u  in  e  n  t 

/*■*/■*  4  '  t-\  t 

U  /  /  V  1  /  OC 

none 

c  a  r  tod  r  aphu 

leaf 

TBD 


Zl-i! 


MG-3P2L/PG 


■  R  0  C  E  3  S 


DESCRIPTION  t 

KEYWORDS  .* 

sources: 

SECURITY  : 
RESP  fd: 
ATTRIBUTE  : 


Thi;.  algorithm  converts  an  i n nut  SecSrsi-iiic  cuordinste 
1  a 1 1 tude/ 1 ons i tude  f'Sir  to  tins  eamvaient  Lo#uei  t/polar  arid 
northing /easting  coordinate  set. 

algorithm 

3.2.115  CGIOBICOA/Part.l 

U 

JUG 

VALUE : 

abbreviation 
ti?e.of.source 
date_acu>Ji  red 
pcocesiiria.done 
m  a  1 1 1  e  m  a  t  i  c  a  i  _  T  i  e  1  d 
tree-level 

I'eccj  i  r  ements-f'e  rf  o  i  niofice 


L  a  i,  /  L  o  1 1  a  _  2  _  L  a  m  U  e  r  t  /  P  o  1  a  r  _  G  r  i  d 

doc  i ih e 1 1 1 
\j  '  /  v  1  /  o 
n  o  n  e 

La  I- tog  I  a  H(i  a 

leaf 

TBD 


2°t- /X 


8  MG-GP2NE 


F  R  0  C  E  3  3 


description: 

keywords: 
sources: 
security: 
resp  pd: 
attribute: 


This,  algorithm  converts  and  inputs  latitudes  lunOi  lode  pair 
into  its  eoui valent  norlliiuS  and  esiiins  Luuruiuale  set* 


algorithm 


3. 2. S3 


CGlOol 00A /Pa rt_I 


U 


JUG 


VALUE 

sbb  rev i at i on 
tyee_of _sou rce 
date.acuuired 
processin^_done 
mathematical-field 
l ree_ 1  eve  1 

r  e  q  u  i  r  e  in  e  1 1 1  s  _  p  e  t-  f  o  r  n  a  r  i  c  e 


Geu  J  r  ipii  i c_2_No  r  th  i riS/Ess  L  i  nsa 
docuee n t 
'07/01/32' 
n  one 

c  a  i- 1  o  a  i  a  p  h  a 

leaf 

TBD 


2  S  w.l 


9  MG-GRDZONECNV 


PROCESS 


description: 

keywords: 
sources: 
security: 
resp  pd: 
attribute: 


This  algorithm  converts  input  longitude  and  latitude 
to  UTM/UP3  Grid  Zone  number  and  letter.  It  fails  at  the 
upper  (closed)  limits  of  both  loriditude  and  latitude. 


alaurithm 

3.2.90  CGIOSIOOA/Par t_I 

U 

JUG 


abbreviation 
tspe.of.sou rce 
d  a  t  e  _  a  c  q  u  i  r  e  d 
processinS-done 
niatheni3tiCel_f  ield 
robus tness 
tree-level 

requirements -performance 
references 


VALUE : 

Grid_Zone_3ene ration 
d  o  c  ij  in  e  n  t 

/  A  ^  »  a  <  / 

'  J  /  /  \J  1  /  O 

a  n  a  1  a  c  e  d /  Pascal  ice  d 
v  a  r  t u a  r  •jrh'j 
3 

lea  f 

Sene  r .  j  1 1  a  _  s  3 1 1 s f actors 
C  JUG  1  CRVGZMG/M6GZD3 


14  MG_L/PG2GP 


P  R  G  C  £ 


description: 

keywords: 
sources: 
security: 
resf  pd: 
attribute: 


This  algorithm  converts  Lawuert/polar  grid  input 
coordinates  into  equivalent  geographic  c  u  u  ruinate 
lati tude/lonsi tude . 


algorithm 

3.2.116  CGiGSlOOA/Part.I 


U 


JWG 


VALUE 

abbreviation 
LaPe_oP_source 
date.acitui  red 
p  recess i ng_done 
mathematical-field 
lree_level 

i' e Ciu  i  r  e in e n t s _ p e  r  f  o  nuance 


Lambe  r t/Po 1  t  r  _3  r  i  d-2 _L a  t  /  Lori 
u’gcu irien  t 
' 07/01/32 ' 
n  o  n  e 

cartography 

1  e  3  r 

TBD 


2  H  -  t  s' 


i 


15  MG-LAMBERTCG 


f 


r R  G  C  E  3 


description: 

This  algorithm  calculates  the  Lambert  constants  i  >.■  mi  ,  e •  J 
for  use  with  the  Lambert /Fulor  Grid  to  UPS  and  the  U F G  to 
Lambert/ Polar  Grid  conversion  algorithms. 


KEYWORDS:  algorithm 
SOURCES:  3.2.11? 

security:  u 
resp  pd:  jug 
attribute: 


1  o  1  0  8 1 0  0  A  /  F'  a  r  t  _  I 


value: 


abb  ro v i at i on 
tsee-of- source 
d  a  t  e  _  a  c  u  u i red 
r  r  o  c  e  s  s  i  n  £  _  d  o  1 1  e 
m  a  t  h  e  m  a  t  i  c  a  1  _  f  i  e  1  d 
robustness 

tree-level 

i  e o u i  re m e n  t  a _ v*  e  r  f  o  r iii a n c >* 
i  e  f  e i enc es 


L  a  m  b  e  r  t  _  C  o  n  s  t  o  n  I  -Gene 
docu me ri  t 
’  0  7  • '  v  1  / '  3  2 

a  n  a  1  s  c  e  d  /  r  a  s  c  si  i  c  d 

c  a  r  t os  r  j .  ■  h  a 

6 

.frd" 

sati  s  f  a  c  t  o  r  -■ 

'  1  JWGjLCjPRO c 


I  c  v  i  u  I  : 


1  cf-  /<, 


I 


1 


17  MG-MAPCUR2NE  PROCESS 

DESCRIPTION : 

Tiiii,  a  1  a  o  p  i  t  i  l  m  converts  an  a  m  a  p  cursor  position  to  its 
eaui valent  northinS/eastinS  ur  1  at i tuue/ 1  una i tuue  or  arid  Lurie 
number /letter  a n d  spheroid  output  coordinate* 


KEYWORDS 
SOURCES 
SECURITY 
R  E  3  P  F  D 


algorithm 

3.2.3? 

U 

JUG 


C 3103100m /Part_I 


attrisute: 


abbreviation 
lype_of_source 
date.a  ccniired 
processinS-done 
m  a  t  h  e  Hi  a  t  i  C  a  1  _  f  ifeld 

tree-level 

re tiui  i  -  enter i t s_pe  r  f  o r nance 


VALUE : 

M  a  p  _  C  u  t  s  u  r  _  2  _  N  j  i 
d  o  c  u  hi  e  n  t 
'07/01/32' 
n  o  n  e 

c  a  r  t  o  3  r  a  p  h  ■< 
iTi  i  d  d  1  e 
TBD 


D 


l  H-  '  7 


1 


10  MG_NE2GP 

description: 


Tiiii-  alSorit'nni  converts  and  inputs  no  r  tii  i  i  id/e-aa  l  *  i  id  set 
lfito  a  1  a t i t u d e /  1  u 1 1 s i 1 1 j i J e  f-air. 


KEYWORDS:  algorithm 
SOURCES:  3.2.34 


C  G 1  0  3 1 0  0  A  /  P  a  c  t  _  I 


security:  u 


R ESP  P D :  J  W G 


attribute: 


value : 


abb  rev  i  at,  i on 
t'ifa.o  f  _sou  rce 
da  l  e  _  a  c  c<  ■  j  x  r  a  <  J 
h  rocessind-done 
a  a  l  n  e  ro  a  1 1  c  a  1 f  i  e  1  d 
l.ree_Ievel 

i  &uui  r e ni e r.  t  s  h  e  r  f  o  s  itsfici 


N  o  i  ■  l  i  i  i  1 1 3  /  E  i  a  t  i  r  < : 
d  o  c  u  in  e  n  t 
'07/01/32 ' 
n  o  n  e 

L  a  !'  t  U  2  I  a  a  h  a 


:  ,.t  ’L 


leaf 

TBD 


MG-NE2MAPCUR 


i 


r r  u  C  E  3  3 


description: 

This*  algorithm  converts  either  no  r  t'ni  ni/SdS  t  i  ns  or 
ialitude/lonaitude  or  -irid  zone  numbe r/let te r  pdii  s  into 
fiinivalent  ;;y  map  cursor  position  pairs  Fur  displas. 


KEYWORDS:  algorithm 
SOURCES.*  3.2.33 
security:  u 


CGlG81vOA/Fert_I 


RESP  F  u :  JUG 


attribute: 

abb  rev i a  1 1 on 
tspe_o  F_sou  rce 
d  a  t  e  _  3  c  u  u  i  r  e  d 
p roeessind-done 
111  3  t  h  e  Hi  a  1 1  C  a  1  _  F  ieid 
tree-level 

r  e  a  u  i  >■  e  in  e  ri  t  s  _  p  e  r  f  o  r  in  a  r  i  c  e 


value: 

No  r  tiling. 'Eas  tiiia_ 

dOCU  ill  e  n  t 

'07/01/32 
ri  o  n  e 

cartography 
ili  i  d  d  1  e 
T  B  D 


description: 

keywords: 

sources: 

SECURITY : 
RESF  fd: 
attribute: 


Th  1  s, 
eoib^osite 


slaorithm  eoriverta  northing 
U T  M  Pair. 


<3 1  id  e  a  a  1 1  i  i  a  iu  a 


alsori thm 

3.2.36  CG10Si00A/p3rt_I 


U 


JUS 


abb  r  ev  i  at  i  or* 
lafe.of.ioijrce 
d  a  t  b  _  a  c  &  u  i  r  e  d 
i--  r  j c e  5  3 1  ri a.done 
in  a  t  h  e  m  3  t  i  c  d  1  _  f  i  e  1  d 
tree. level 

r  e  q  u  i  r  e  ih  e  1 1 1  s  _  ^  e  r  f  o  r  iu  «  n  c  e 


VALUE : 

f  <  Cl  I  o  I  n  1 1  Si  /  E  o  5  t  X  ;  ■  a 

docu Hi e i'i  t 

v  /  /  'j  1/  o*; 
ri  One 

Cal  tea  I  a  I-  ha 
lea  r 
T  B  D 


2*1-1  o 


1 


r  r.'  j  C  E  3  * 


i 


A 


3  MG-PG2UPS 


FRGCE3S 


description: 


keywords: 
sources: 
security: 
resp  pd: 

ATTRIBUTE  : 


This  alsorithia  converts  F'uldr  Grid  no  r  t  i  ■  i  ni/sds  t  i  nd 
coordinates  into  eauivsient  Universal  Polar  3  te  i  os  rarh  i  l 
ur  Universal  Transverse  Mercator  (utilicinS  the  NE2UTM 
algorithm)  . 

alsorithiii 

3.2.118  CGI  0810 Oh/ Part-1 

U 

JUG 

abbreviation 
l  s  p  e  _  o  f  _  s  o  >J  r  c  e 
Uate_scQui red 
p  r  o e  e  s  s  l  n  s  _  do n e 

Hi  a  1 1 1  e  hi  a  l  i  C  a  1  _  iA  i  e  1  d 

robustness 
tree-level 

r  e  ct  u  i  re  in  e  n  t  s  _ e  r  f  o  r  m  a  i'i  c  a 
ref e r  e  1 1  c  e  s 


value: 

P  o  1  a  >■  -  G  r  1  d  _  2  _  u  F'  3 
d  o  c  ij  in  e  n  t 

'  0  ?  01/32  ' 

an elated /Pascal  iced 

a  i  t  u  a  r  i  p  i'i  a 

leaf 

TUB 

il  J  W  u  3  F  3  2  U  F  3 


| 


2  9-21 


27  MG-RAD2DE3 


PROCESS 


DESCRIPTION*. 


This  algorithm  converts  an  ansle  in  scaled  Pi  radians 
into  an  eauivalenu  ansle  iri  u»;S  fees  /  m  irm  tes/seconds  « 


SOURCES*.  CG108100A/Part_I 


security:  u 


resp  pd:  jug 


attribute: 


value: 


abb  rev  1  at  1  on 
tape_of_source 
i J  a  t  e  -  a  c  a  u  i  r  e  d 
processing -do n  e 
matheniatical_f 
tree-level 

i-  e  u ij  t  r  e ni e ri  t s  _ p  e  r  f  o  r  m  a ri  c  e 


Radians_2_DeSiee-. 
docy  in  e  n  t 
'  0  7  /  0  i  /  3  2  ' 
c  i  o  n  e 

t  T  1  a  0 1 1  o  ill  e  1 1' 

lea  f 


Z  9-2  Z. 


3  MG_UPS2PG 


description: 

keywords: 

sources: 

SECURITY  : 
RESP  fd: 
attribute: 


PROCESS 


This  alaorithm  converts  Universal  Polar  3  Leroa  rapine  mlo 
eouivalent  polar  arid  coordinates* 


a  1  a  o  r  i  t  h  ru 

3.2.117  CG10S1C0A /'Pa  rt_I 

U 

JWG 


abb  rev i at i on 
tspe_o  f_  sou  rce 
date_acc<uired 
p  rocessi n a -done 
mathematical-field 
t  ree_l e ve 1 

rectui  rements-performance 


value: 

U  P  3  _  2  _  P  u  1  a  r  _  G  r  i  >  j 
d  o  c  u  hi  e  n  t 

/  A  /  A  <  <  o  n  / 

Mi/  Ml  /  d  *- 

n  u  n  e 

C3 I toai aphu 
leaf 
T  B  D 


34  MG-UTM2NE 


PROCESS 


description: 

keywords; 
sources: 
security: 
resp  pd: 
attribute: 


This  algorithm  converts  a  UTM  coaidumU  set  into  tiie 
equivalent  composite  nor  tiling  and  e.istinS  r.ji  r , 


algorithm 

3.2.85  CGI vSlOOA/Part.I 

U 

JWG 


abb  rev i at i on 
Lsfe.of.source 
d  a  t  e  _  a  c  c*  u  i red 
k rocess i nS-done 
mathematical-field 
t  ree-level 

r  e  « i  j  i  r  e  m  e  n  t s  _  p  e  r  f  o  r  m  a  n  c  e 


VALUE : 

UTM -2-No  r  th  i  nS/Ea  s  t  i  ni 
d  o  c  u  m  e  n  t 
‘  0  7 -*0 1/32  ' 
n  o  n  e 

c  a  r  t  o  5  r  a  f-  h  a 
leaf 


6400 

6500 

6600 

6700 

6800 

6900 

7000 

7100 

7200 

7300 

7400 

7500 

7600 

7700 

7800 

7900 

8000 

8100 

8200 

8300 

8400 

3500 

8600 

8700 

8800 

8900 

9000 

9100 

9200 

9300 

9400 

9500 

9600 

9700 

9800 

9900 

10000 

10100 

10200 

10300 

10400 

10500 

10600 

10700 

10800 

10900 

11000 

11100 

11200 

11300 

11400 

11500 

11600 

11700 

11800 

11900 

12000 

12100 

12200 

12300 

12400 

12500 

12600 
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URITELN  <'  it iu 


INPUT  PARAMETERS}- 


0.25; 

0  ♦  0  J 
0.125J 

o.o; 

0  ♦  0675 ; 
0.375; 


setup  Em-fiiir  j  ui  iu  j  ■ 
Si.’iui*  Sum  i  mu-ju  r  Art  L  a 


icai'iinjL'n’j/i 
t Sum i mu Ju i  An  La)  ? 


URITELN  ( 
URITELN  < 
URITELN  < 


S e  t u P  LumuuuL  1  a 
auiun  LumbuoR  La 
at-  tup  P'n  i  U  in  r 
au  iup  Pii  L  L  i  a  '  r 
it'  iur  P  i  i  i  1  i  a  r 
au  iu;-  Pi  i  L  2  i  a  t 


a  '  r  L  tiuLniil.  )  ; 
a  '  >  I.  tiwi'J .1 R  )  J 
r  Pi i  i u ) ; 
tPiiiL) ; 

»Piiii  > ; 
t  F  i  i  L  2  )  ; 


L -n»ibn  r  iCm  ia  i<in  iGt.'iu.'  i' <i  i  i  m  i  ( 


■C  lit- 1  uuir  ui  r  ui  c-eettra  f  r 

J  C  -  X 


■CGI.0BAI  } 

{GLOBAL} 

<GUT> 

{OUT} 

{OUT} 

{OUT} 

<GUT> 

<DUT> 

{IN} 

{IN} 

<  I  N> 

{IN} 

■C  TN> 

•C I  N> 

urn  lt.-ii} 


F  m;  ..'II  i  i.  V  i  i  a  * 

St* hi  i  Miuviu  rfix  i  a  > 

LumbuuL * 

Lemb'ioR  > 

Pi i  l  U  t 
PiiiL » 

F  i  i  i  1  t 
Piii2  > 

K  d  ? 

Iul.il 

P  r uJtL  iiunCuiiuRduiui  i 
LembubC  f 

Suuu  i  hfiiEi.'i.  u'ii  i  i  i L  i  a  t 
He* i ai  iie  i  t  )  J 


ri  rS  rS  rS  r  ■>  rS  fS  /S 


12700 

•C 

s 

12800 

WRITELNi 

12900 

WRITEI.H 

\  ' 

K-ii-’i'j  i  j  '  >  K  j;-';-'  i  )  t 

13000 

WRITELN 

(  ' 

I  L  <J  '  1  lulii)  ! 

13100 

WRITELN 

<  ' 

P  I'U  JliL  L  iui  iC'Jlll-T  R  JI.I  L  U\>  i.  i  *  7  P  I'J  L  i  unClJI  H.'R^U  i  IJj  >  » 

13200 

WRITELN 

(  ' 

l.ciUlU'jclC  A  ii  '  f  L  )  » 

13300 

WRITELN 

<  ' 

SullOI'li'jctLtfull-  Lc  l  l':<  i'j  '  7  3  U  ■  J  -i  f  Si*  JE  U  ■?  II  L  l‘  i  C  l  L  "1  7  7 

13400 

WRITELN 

<  ' 

Htniiui'iifie  j i  Ht.’ in  i  !»f  i  if  i'U' ) 

13500 

-c 

** 

J 

13600 

END. 

■C  D  r  i  ve  rLsiHue  i 

r  LCuiii>  Lt;ii  LGei it  res  L  Ion  PROCEDURE 

100 
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w  O  ‘  J  •J 

2  6  0  0 
2  2  0  0 
2  8  0  0 

2  7  0  0 

3  0  0  0 
3100 
3200 
330  0 
3  400 
3500 
3600 
3700 


PROCEDURE  Conve  r  tToUF'S  p  e  r  f  o  r  iti  s  the  P  o  1  a  r  Grid  to  Universal 
Polar  3 1  e  r  i  os  r  aph  i  c  conversion  algorithm  described  in 
paragraph  3.2.118  of  CG108100A  dated  23  October  1978. 

PROCEDURE  ConvertToUF'S  is  CALLed  by!  F'ROCEDURE  F'ularTuUr'3 

PROCEDURE  F'olarToUTM 


F'  r  o  d  r  a  hi  Hi  e  d  by  J.w.dillis 


5-6-32 


This  procedure  ASSUMES  that  certain  data  are  available  as 
required  by  the  algorithm  but  nut  au'eauately  described  in  tile 
algorithm  description. 

This  procedure  DOE3NOT  perform  data  validity  checks  that  are 
not  specific d  in  the  algorithm  description.  This  is  to  allow 
the  algorithm  features  to  be  presenters  mure  clearly. 


3700 
•1000 
1100 
1200 
1300 
1-100 
1500 
1600 
1700 
1300 
170  0 
5000 
5100 
5200 
5300 
5100 
5500 
5600 
5700 
5300 
5700 
6000 
6100 
6200 
6300 


F'ROCEDURE  ConvertToUF'S  accepts  polar  drid  northind  and  easting! 
coordinates  and  con  verts  them  to  UF'3  coordinates.  I 


Meters 

Decameters 

Letter 


I  n  d  e 

3  r i dL  ei te  r 


-  REAL? 

=  REAL? 

-  '  A  '  .  .  '  Z  '  J 


j  integer; 

;  ARRAY  Cl.  .263  OF  Letter 


FUNCTION  a M 0 D b  (<IN>  a , b J REAL ) t INTEGER J 

r 

VAR  Ainteser  1  INTEGER? 

Binteser  J  INTEGER? 

•C 

BEGIN 

AiriteaerJ-=TRUNC(a)  ? 
Binteaer5=TRUNC(b) ? 
a  M  0  D  b  J  =Ainteder  MOD  Binteder 
END?  Oof  FUNCTION  M0Dab> 

0 

< 

J  O  -to 


i 


rS  f'-\  rS  rS  cS  rS  *"S  rS  rS  rS 


s 


6  40  0 
6500 
66  00 

6  700 
6800 
6900 
7000 
7100 
7200 

7  300 
7400 


BEGIN 


{Initialize  the  Grid Letter  array.  > 

{  i 

FOR  Inde ; ;  1-  1  TO  26  DO 

G  r  i  dLe  tte  r  C  I  nde;-;  1 1  =  CHR  (  ORD  (  A  '  )  + 1  n d e —  1 )  5 
{  > 

{  7 

{  > 

{Calculate  the  UF’SEast  i  naN'jmbe  r-  %  UF'SEa st  i  nsLet  te  r  I  ride;;  7 


IF  F'GZuneLette  r  >  'M' 

THEN  BEGIN 

U  F  S  E  a  s  t  i  n  a  N  u  hi  b  e  r  t  —  a  H  G  D  b  (  (FGEsstinaCoord-ZOOOOO)  . 

100000)/ 10 ; 

I  ride;;  5  -=aM0Db  (  <  FGEas  t  i  n^Coo  rd-200000  )  / 1000000  .  20  ) 
END!  {of  I F  7 


9000 
9100 
9200 
9  30  0 
9400 
9500 
9600 
9700 
9800 
9900 
10000 
10100 
10200 
10  300 
1 0  400 
10500 
106  0  0 


IF  P G Z u n e L e 1 1 e r  <  ' N  ' 

THEN  BEGIN 

UFSEas  t  i  naNuiribe  r  !  -  aMODb  (  FGEas  t  i  n3 Coo  r  d  » 

100000)/ 10! 

I  n d e ;; !  =aM0Iib(  FGEas  tineiCoord/1 00000  » 13) 
END!  {of  IF7 

r 

UR  I  TEL N  ? 

URITELN  (  ' UPSEASTINGNUMBER  IS  '  .  UFSE AST  I NGNUMBER  )  ! 
URITELN  ('AND  ITS  LETTER  INDEX  IS  '.INDEX)! 

{ 

{Calculate  the  UF'SEastingLetter  from  its  index. 

{ 

IF  Inde;-;  >  16 

THEN  Inde;;  J  = Index+2 ! 

IF  (Inde;;  >  2)  AND  (Inde;;  <  17) 

THEN  Index  J - Index+ 1 ! 

UFSEas  l  inaLette  r  5  =  Gr  idLette  t-C  IndexI ! 

{ 

URITELN ! 

URITELN  (  ' UPSEA3T INGLETTER  IS  ' » UPSE AST  I NGLETTER ) ! 
URITELN  ('AND  ITS  CORRECTED  INDEX  IS  '.INDEX)! 


1  0800 
10900 
11000 
1 1 100 
11200 
11300 
11400 
11500 
11600 
11700 
11300 
11900 
12000 
12100 
12200 
12300 
12400 
12500 
12600 


{Calculate  the  UFSNo  rthirnSNumbe  r  S  Nu  r  tii  i  na'Le  t  te  r  Index  > 

{  1 
IF  F’GZoneLetter  >  'M' 

THEN  BEGIN 

U  P  3  N  o  r  t  h  i  n  si  N  u  m  b  e  r  J  -  a  H  0  B  b  (  (  P G No  r  t h  i  n U C u o  r  d  -  i  3 0 0 0  0  0 

i 00000 ) / 1 0 ! 

I  n  d  e  1  —  a  M  0  D  b  (  (  (  F'  G  N  o  r  t  h  i  n  a  C  o  o  r  d  - 1 3  0  0  0  0  0  )  /  1  0  0  0  0  0  )  ? 
24) 

END!  {of  I F> 

r  •> 

A.  S 

IF  FGZoneLetter  <  'N' 

THEN  BEGIN 

UFSNo  r  th  i  naNumbe  r  I  =aM0Db  (  (  F G N o  r  t i  i  i  n 5 C o o  r  d -  3 0 0 0 0 0  ) 

100000) /IQ! 

I  n  d  e ;;  5  -  a  M  0  D  b  (  (  (PG  N  or  thindCuord -  800000)/ 10  00005  > 
24) 

END!  {of  IF> 

{  7 

{  1 

lo-;/ 


12700  WRITELN; 

12S00  WRITELN  (  ' UFSNGRTH I NGNUMBER  IS  '  » UP3NQRTHINGNUMBER )  i 
12900  WRITELN  ('AND  ITS  LETTER  INDEX  IS  ', INDEX); 

1.3000  -CCe  1  cu  1  3  t  e  the  UPSNo  r  til  i  ndLe  1 1  e  r  from  its  Inde;; 

13100  -C 

13200  IF  I  n  d  e ;;  1  6 

13300  THEN  I ncie;;  •  =  I  ricie;;  +  2  > 

13400  IF  Inde,:  <  17 

13500  THEN  INDEX  :  =  I r.de;;+  1  » 

13600  U  F’ SNort  in  iridLetterl^Grid  Letter  Elude;;  3) 

137  00  \ 

13300  WRITELN) 

13900  WRITELN  (  '  UPSNORTH I NGLETTER  IS  '  >  UF'SNORTHINGLETTER  >  ; 
14000  WRITELN  ('AND  ITS  CORRECTED  INDEX  IS  '» INDEX); 

14100  < 

14200  < 

.14300  \ 

14400  < 

14500  END?  Oof  PROCEDURE  Conve r tToUPS> 

14600  < 

14700  END.  Oof  MODULE  CNV2UPS> 


3  O  -  t  X 


TEMPORARY 


100 

200 

300 

400 

500 

600 

700 

aoo 

900 

1000 

1100 

1200 

1300 

1400 

1500 

1600 

1700 

1300 

1900 

2000 

2100 

2200 

2300 

2400 

2500 

2600 

2700 

2300 

2900 

3000 

3100 

3200 

3300 

3400 

3500 

3600 

3700 

3800 

3900 

4000 

4100 

4200 

4250 

4900 

5000 

5100 

5200 

5300 

5400 

5500 

5600 

5700 

5800 

5900 

6000 

6002 

6100 

6200 

6300 

6400 

6500 


STUB  PROCEDURE 


MODULE  CNV2UTM  < I NFUT » OUTPUT > » 

< 

< 

TYPE 

Meters  =  RE 

Decameters  =  RE 

Letter  =SE 


apiieres 


=REAL; 

=real; 

=  SET  OF  CHAR, 
= INTEGER , 


PROCEDURE  ConvertToUTM 


(•CIN  >  PGNorthingCoord 

■CIN  >  FGEastinSCoord 

•CIN  >  FGZoneLette r 

CIN  >  F'GZoneNumbe  r 

CIN  >  Spins roidNumber 

C0UT>  VAR  UTMZoneLette  r 
CQUT>  VAR  UTMEastinsLette r 
C0UT>  VAR  UTMNorthinSLetter 
C0UT>  VAR  UTMEastinSN umber 
C0UT>  VAR  UTMNo r th l nSNumbe r 


t  Mete  r s  5 
: Me  te  rs  f 
1  L  e  1 1  e  r  • 

•.integer; 

t  Spheres  f 
;  L  e  t  te  r  5 
•Letter; 
JLetterJ 
;  Decameters; 

; Decameters  - ? 


C  PROCEDURE  ConvertToUTM  performs  the  Polar  Grid  to  Universal 
C  Transverse  Mercator  conversion  algorithm  described  in 
C  paragraph  3.2.86  of  CG103100A  dated  23  October  1978, 


C  PR 

C 

C 

C  Pr 
C 

C  Th 
C  re 
C  al 
C 

C  Th 
C  no 
C  th 
C 

C  PR 
C  co 
C 

c 

TYPE 


PROCEDURE  ConvertToUTM  is  CALLed  bat  PROCEDURE  Polar To UTM 

PROCEDURE  PularToUPS 


Proa  rammed  bs  J.w.gillis 


5-6-82 


This  procedure  ASSUMES  that  certain  data  are  avails ble  as 
reeuired  by  the  algorithm  but  not  adeouatels  described  in  the 
algorithm  description. 

This  procedure  D0ESN0T  perform  data  validity  checks  that  are 
not  specified  in  the  algorithm  description .  This  is  to  allow 
the  algorithm  features  to  be  Presen tes  more  clearly. 

PROCEDURE  ConvertToUTM  accepts  polar  grid  northing  and  easting 
coordinates  and  converts  them  to  UTM  coordinates. 


Meters 

Decameters 

Letter 

Spheres 


*  real; 

=  real; 

=  set  of  char; 
=  integer; 


BEGIN 

<  IF  > 

■CTHEN> 

URITELN  ('  CNV2UTM  HOOKED  UP  0K'> 

< 

< 

END?  lof  PROCEDURE  Cor.ve  r  tToUTM> 

< 

END.  <of  MODULE  Cnv2UTM>  Jo-J-i 


*i  rS  rS  rS  f'i  rS  r~>  r'*i  r‘*»  r -I  r-%  r\  rS  t* i  r*>  r'n 


too 

TO 

300 
•4  0  0 

PROGRAM  ConverLToUTM  (I 

to 

TYPE 

N  r  U  T .OUTPUT) ? 

5  0  0 

ASCIIArraa  -A 

RR AY  Cl. *41  OF 

CHAR  7 

600 

700 

300 

Letter  =  CHAR} 

r  ^ 

~y.j 

V  A  R 

n  /•>  r\ 

7  V  '.,*' 

UTMCoo  rd i na  te 

t  RECORD 

1000 

U  T  M  Z  o  n  e  N  u  in  b  e  r 

t  ARRAY C 1 . .23 

uF  INTEGER? 

t  too 

UTMZoneLe  tte  r 

t  L  e  1 1  e  r  5 

1.200 

UTMEastindLetter 

t  L  e  1 1  e  r  5 

.1300 

UTMNorthinSLetter 

•Letter? 

t  400 

UTMEsstiridNu  m  b  e  r 

t  A  R  R  A  Y  C 1  .  .43 

OF  CHAR  5 

1500 

1300 

U  T  M  N  o  i  t  i  i  i  n  S  N  u  hi  b  e  r 
END?  -CUTMC  0 

: array: i. .43 

OF  C H A R 

1 3  0  0 

F'GNorthingCoord 

t INTEGER? 

T I N  0 

i  ?  0  0 

P G  E  a  s  t i n  gCuor  d 

t INTEGER? 

TIN  0 

2000 

FGZoneLetter 

(Letter  ? 

-CIN0 

2  i  0  0 

PGZoneNufiibe  r 

t INTEGER? 

TING 

220  0 

3  p  h  e  r  0 i d N u  m ber 

t INTEGER  ? 

TING 

2  2  0  0 

I*  EletterOf f set 

,  *  Tel*  T  p f 

Tnl  t 

2  -t  0  0 

M  u  iti  b  e  r 

*.  INTEGER 

J 

2500 
2  6  0  0 

2300 

2  V  0  0 

3  0  0  0 
3  tOO 
3300 
3  3  00 
34  00 
3  5  0  0 
3  6  0  0 

3  7  0  0 
3300 
7>  ?  0  0 

•  *  a 
••  •  •; 

"  *.  !J  0 

■4  7  0  0 
430  0 

4  4  00 
■4  500 
■4300 
4700 
■4  300 

4  V  0  0 

5  C  0  0 
5100 
5  2  0  0 
5  3  0  0 
5-400 
550  0 
5300 
5  7  00 
5300 

5  7  7  0 
0  0  0 

it  00 

6  200 
'300 


T  P  R  0  C  E  D  U  R  E  C  o  n  v  e  r  t  T  o  U  T  M  pii-forms  the  r  u  i  a  r  Grid  to  Universal 
T  Transverse  Mercator  conversion  algorithm  described  in 


T  PROCEDURE  CorivertToUTM  is  C  A  L  L  e  d  bat  PROCEDURE  F'oiarTuUTM 
-C  r  R  0  C  E  D  U  E  E  P  o  1  a  t  ■  T  u  U  P  3 


-C  This  procedure  ASSUMES  that  certain  data  are  available  as  > 

<  reouirud  ba  the  algorithm  but  not  adequately  described  in  the  3- 

{algorithm  description*  3- 

\  F 

T  This  procedure  D0E3N0T  perform  data  validity  checks  that  are  > 

-C  riot  specified  in  the  algorithm  description.  This  is  to  allow  G- 
-C  the  3 1  0  o  r  1 1  h  iii  features  to  be  presences  more  c  1  e  a  r  1  a  .  > 

TO 

T  PROCEDURE  ConvertToUTM  accepts  polar  arid  northing  and  eastingO 

<  coordinates  and  converts  them  to  UTM  coordinates*  3 

TG 

P  R  0  C  E  D  U  F;  E  C  o  n  v  e  r  t  T  o  A  S  C 1 1  (  Number*  INTEGER  T  i  n  0  ? 

VAR  A 1 r haNum « A3C 1 1  A r r as  5  5 
r 

J 

T  Convert  an  integer  number  to  its  ASCII  rep  re  sen  ta  t  ion  in  base  1 
T  This  Procedure  not  detailed  in  source  0 

■c  o 

V  A  R  It  INTEGER? 

A  d  J  t  INTEGER? 

BEGIN 

PGR  I  t=  4  D0WNT0  1  DO 
BEGIN 

idj  5  -  (Number  MOD  10)  t  G R D ( ' 0 ' 5  ? 

A  1  p  h  a  N  u  m  C  I  3  S  -  C  H  R  (  A  d  J  /  ? 

Number  t =  Number  DIV  10 

END 


3  &  ■ '  h 


C  C  o  ri  v  e  r  t  T  o  A  3  C 1 1  > 


BEGIN  -CC onve  r  tT oUTM> 

-c  y 

x  >'  i?  d  o  i.  i  i  i '  ■  j  C  s 

-c  y 

writeln; 

WRITELN  i.  '  F Q  Northing  CuOi'd  (integer)  > 

R  E  A  Ei  L  N  (  P  G  M  o  r  L  h  i  1 1  d  C  u  o  r  d )  5 

WRITELN  i  'FG  Ess  tins  Coord  (integer)  ')  ? 

READLN  (  F"  G  E  a  s  L  i  r  i  -i  C  o  u  r  d )  5 

WRITE  L  N  <  '  F'G  Zorie  Letter  (Letter)  )  5 

READLN  <  F’GZoneLet  Le  r  )  5 

WRITELN  ('F'G  Zone  Number  (inteder)  ')  t 
READLN  (PGZoneN umber)  5 

WRITELN  <  '  3  p  h  e  r  o  i  d  N  u  iiiber  (  i  n  teaer)  '  )  ? 
READLN  (  Sr-he  r o  i  dNumbe  r  )  * 

r 

S 

x  end  of  input  > 

f  T 

\  s 

WITH  WTMCoordinate  DO 
BEGIN 

•I  Find  Eastind  Letter  (A2)  > 

I  1  =  p  G  Z  u  n  e  N  u  in  b  e  r  M  G  D  3  * 

■C  Eastind  Letter  Offset  Factor  > 

CASE  I  OF 

0:  ELetterOf fset  5-  24? 
i:  ELetterOf fset  1  =  65 

2?  ELetterOf  fset  *.  =  15 
END? 


?  6  0  0 
9  70  0 
9  E  0  0 
9900 
10000 
10100 
1  0200 
103  0  0 
1 0  4  O  o 

10500 
1 0  6  0  0 
10  700 
10800 


Tel  5-  (  PGEa s t i nd Coo r d  -  500000  )  DIv  100000  +  ELe t te ru f fae t ? 
•CO 

-C  Following  letter  conversion  not  according  tu  duwumen  ta  t  i  on  ?  > 
x  Procedure  in  documentation  fails  at  l«sst  t  8®d  Hersfeld*  FRG» 

-C  which  is  in  32U>  and  this  one  succeeds  at  least  there.  > 

CO¬ 
IF  Tel  13  THEN  Tel  1  =  Tel  +  2 

ELSE  IF  Tel  >  10  THEN  Tel  1  =  Tel  -  1 
ELSE  Tel  » -  Tel  -  2  ? 

Tel  t  -  Tel  +  ORD('A')  -  1  ? 

UTMEastindLetter  1  -  C  h  r  <  Tel  )  -C  A  2  >  ? 

-C  F  i  n  d  Northi n d  Number  > 

IF  NOT  0 D D (  P G Z o n e N u m b e r  )  THEN 

F  GNo  r  til  i  ndCoo  rd  1  ~  F' G N o  r  t h  i  n d C o o  r d  +  50  000 0 » 


1  1  L  0  0 

T  P 

PGNo  r  th iudCoo i d  MOD 

i.  1.200 

»  ft 

1  : =  Tp  DIV  100000  +  1? 

1 1  3  0  0 

\ 

Hake  Spheroid  Adjustment 

1 1  400 

CASE  3 p h e r o i d N u m b e r  OF 

1 1.300 
116  00 

•C  Clark  1366  > 

1 1 "  0  0 

11  IF  (  F‘ G Z o n e N u m b e  r <  3 1 

11300 

ELSE  IF  ( POZoueNijiKue 

1.1900 

THEN  T ci  1 

12000 

<  International  > 

12100 

2  5  IF  ( PGZoneNUMBER  1 46 

12200 

12300 

c  y 

THEN  Tri  1  l  ~  Tnl  + 

1 2  4  00 

■C  C I 3  r k  1330  > 

12500 

12600 

/*  ■» 
\  ■* 

31  Tnl  1 r-  Tnl  +  10“ 

Q  !  intll  I  I  I  A  .  -  i  n  1  T  l 
i  33)  Ctupo  in  soul  cel 


3o-/-r 


i. y  o  o 

12300 
1290  0 
I  3  0  0  0 

1  v-".'  x  y  v 

13200 
13300 
13400 
13500 
13600 
13700 
13300 
13900 
1  4000 
14100 
14200 
14300 
14  400 
1 1500 
1  4  60  0 
1  4  700 
1  4800 
14  90  0 
1  5  00  0 


■C  Everest  > 

4  1  IF  F G Z o i i  e N u iii b e  r 


46  THEN  Tnl  :  =  Tr.l  r  1> 


mth  i  1 1  X 


o 


•C  Bessel  > 

51  IF  (  FQZo i "i eNu m b e  r < 52/  AND  <  FGZoneLe  t  te  r> '  R 
•C  Australian  > 

6  5  IF  0 D D  (  F'  6  Z  o  n  e  N  U  m  b  e  r  )  THEN  Tnl  5  -  T n  1  15 

ELSE  Tnl  1=  Tnl  +  5 

<> 

END  >  0  CASE  G  F  S  p  h  e  r  o  i  d  N  u  iri  b  e  r  3- 

r  ^ 

\  jr 

C  Calculate  Final  N  o  r  t  h  i  n  g  Letter  I  n  d  t- ;;  > 

O 

IF  Tnl  >  14  THEN  Tnl  5  =  Tnl  +  2 

ELSE  IF  Tnl  >  9  THEN  Tnl  5 =  Tnl  r  1? 

Tnl  1 =  T  n i  +  0  R  D ( ' A  *  )  -  1  * 

UTMNorthinsiLetter  1=  CKR(  Tnl  >* 

\  F  o  r  ru  at  for  Gut  ?  u  t  > 

Co  nve  r  t  To  A3C  1 1  (  F'GNo  r  til  i  n  si  Coo  rd  MOD  100000?  UTMNor thing Numbe r  >* 
Con vertTo ASCII!  F'GEast  i  niCuo  i  d  MOD  100000*  UTHEas  t  i  nsiNuuibe  r  )* 
UTMZoneNijBibe  r  C  1  j  5-  FGZoneNuiribe  r  D I V  10* 

UTMZoi  leWijnibe  r  C  2 1  1=  FGZoneNumbe  r  MOD  10  * 

UTMZoneLetter  1-  F'GZoneLetter  * 


X  1  'j  y 

o: 

1 3  2  0  0 

<  WRITE 

OUTPUT 

■% 

/ 

15300 

*•  •% 

\ 

154  00 

WRITELN 

r 

15500 

UR  I  T E L  N 

t 

15600 

WRITELN* 

157  0  0 

.  FOR  I  := 

1  TG 

2  DO  WRITELN  (' 

*  i—  m  ^ 

X  J  O  ‘J  \J 

WRITELN 

(  'UTM 

Zone  Letter 

t  Z  7  0  0 

WRITELN 

<  'UTM 

Easting  Letter 

16  000 

WRITELN 

(  'UTM 

Northing  Letter 

16100 

FOR  I 

1  TO 

4  DO  BEGIN 

loZOO 

WRITELN 

<  'UTM 

E  a  s  t  i  n  si  Number 

16300 

WRITELN 

(  'UTM 

Northing  Number 

16400 

END  ; 

16500 

r 

j" 

1  6600 

<  END  OF 

0UTFU 

T  > 

16700 

v  y 

I  6300 

END  -Cuf 

WITH  U 

TM  > 

16900 

V  } 

l  7  0  o  0 

END.  -CConve  r  tTuUTM> 

U  T  M  Z  u  n  e  N  u  ih  b  e  r  '  *  I  *  U  T  M  Z  o  n  e  N  ■  j  ii<  b  e  i  u  I  j  / 
' * UTMZoneL _ tte r )  * 

*uTM East! n ^Letter)  * 

'  *UTMNorthinsi  Letter*  j 

*  I >  *  UTMEas  t i ngNumbe  r  C  I  j  )  * 

»  I  *  '  *  U  T  M  N  o  r  t  h  i  1 1  si  N  u  m  b  e  r  l  I  j  /  * 


PROCEDURE  ConvertToUTM  IS 


TYPE  ShortArray  IS  ARRAYC1..4)  OF  CHARACTER; 

PGNorth ingCoord 
PGEast ingCoord 
PG2oneLetter 
PG2oneNumber 
Sphero i dNumber 

UTM2oneNumber  :ARRAY<1..2)  OF  INTEGER; 

UTM2oneLetter  : CHARACTER; 

UTMEast ingLetter  : CHARACTER; 

UTMNorth ingLetter  .‘CHARACTER; 

UTMEast ingNumber  :ARRAY(1..4)  OF  CHARACTER; 

UTMNorth ingNumber  :ARRAY(1..4)  OF  CHARACTER; 

1/  EletterOf fset,  Tel/  T p,  tnl.  ASCIIPos:  INTEGER;  — Local 

PROCEDURE  ConvertToASCII (  Number:  in  INTEGER; 

AlphaNum:  out  ShortArray)  IS 

—  Convert  an  integer  number  to  its  ASCII  representation  in  base  10 

—  This  procedure  not  detailed  in  source 
BEGIN 

FOR  I  in  reverse  1  ..  4  LOOP 

ASCIIPos  :=  Number  MOD  10  +  CHARACTER ' POS (' 0 ') ; 

AipnaNum(i)  :=  character'Ualc  asciipos  ); 

Number  :=  Number  /  10; 

END  loop; 

END  ConvertToASCI I ; 

FUNCTION  ODD (  I:  INTEGER  )  return  BOOLEAN  is 
BEGIN 

IF  I  =  2  *  <  1/2  )  THEN  return  false; 

Else  return  true; 

END  if; 

END  odd; 

BEGIN 

--  Find  Easting  Letter  (A2) 

I  is  PG2oneNumber  MOD  3; 

—  Easting  Letter  Offet  Factor 
CASE  I  IS 

WHEN  0s  >ELetterOffet  :s  ?4; 
when  is >ELetterOf fet  :s  6; 

WHEN  2s >ELetterOf f set  :s  is; 

WHEN  OTHERS  s>  NULL; 

END  case; 

Tel  :s  (  PGEast ingCoord  -  500000  )  /  100000  +  ELetterOf f set ; 

IF  Tel  >  14  THEN  Tel  :=  Tel  +  2; 

ELSIF  Tel  >9  THEN  Tel  :s  Tel  +  l; 

END  IF; 

Tel  :s  Tel  ♦  CHARACTER 'POS< ' a '  )  -  l; 

UTMEast ingLetter  : s  character 'UAL (  Tel  >;  — A 2 

--  Find  Nortning  Number 

IF  NOT  OddC  PG2oneNumber  )  THEN 

PGNorth ingCoord  :s  PGNorth ingCoord  +  500000; 

END  IF; 


: INTEGER;  — IN 
: INTEGER;  --IN 
.‘CHARACTER;  — IN 
: INTEGER;  — IN 
: INTEGER;  — IN 


Tp  :=  PGNorth ingCoord  MOD  2000000; 

T n 1  : =  Tp  ✓  100000  +  l; 

—  Make  Spheroid  Adjustment 
CASE  Sphero i dNumOer  IS 
—  Clark  1866 

WHEN  1  =  >IF  PGGr idzoneNum  <31  OR  PGZoneNumoer >  58  THEN  Tnl  := 
ELSIF  PGZoneNumoer  >  51  AND  PGZoneNumber  <  59  --typo 
THEN  Tnl  :=  Tnl  -  10; 

END  IF; 

—  International 

WHEN  2=>IF  PGZoneNumber  >  46  AND  PGZoneNumber  <  52  THEN  Tnl 
END  IF; 

—  Clark  1880 

WHEN  3=  >Tn 1  :=  Tnl  +  10; 

—  Everest 

WHEN  4=>IF  PGZoneNumber  <  46  THEN  Tnl  :=  Tnl  +  10;  END  IF; 

—  Bessel 

WHEN  5=>IF  PGZoneNumber  <52  AND  PGZoneLetter  <'R'  THEN  Tnl 
END  IF; 

—  Austra 1 i an 

WHEN  6=>IF  ODD <  PGZoneNumber  )  THEN  Tnl  :=  Tnl  +  is; 

END  IF; 

WHEN  OTHERS  s>  NULL;  —  Not  in  source 
END  case; 

—  Calculate  Final  Northing  Letter  index 
If  Tnl  >  14  Then  Tnl  :=  Tnl  +  2; 

ELSIF  Tnl  >9  THEN  Tnl  :=  Tnl  +  i; 

END  if; 

Tnl  ::  Tnl  +  CHARACTER 'POSC 'A' )  -  l ; 

UTMNorth ing Letter  ::  CHARACTER 'UAL <  Tnl  )! 

--  format  for  output 

ConvertTo ASCII (  PGNorth ingCoord  MOD  100000*  UTMNor th i ng Number  ); 
ConvertToASCII (  PGEast ingCoord  MOD  100000*  UTMEast ingNumber  ); 
UTMZoneNum ( 1 )  :=  PGZoneNumber  ✓  10; 

UTMZoneNum<2)  : =  PGZoneNumber  MOD  10; 

UTMZoneLetter  :=  PGZoneeLetter ; 

END  ConvertToUTM; 


Tnl  +  10; 
in  source 


=  Tnl  +  10 


=  Tnl  +  10 


i  o  -  '  y 


100 
200 
300 
400 
'300 
A  00 
700 
800 
900 

iooo 
1100 
120G 
1300 
1400 
1500 
1600 
1700 
1300 
1900 
2  0  0  0 
2100 
2200 
2300 
2  400 
2500 
2600 
2700 
2300 
2900 
3000 
3100 
3200 
3300 
3400 
3500 
3600 
3700 
3800 
3900 
4000 


1 

PROGRAM  DRVGZTB  (  INPUT , OUTPUT  )? 


•1  PROGRAM  DRVGZTB  provides  <3  test  driver  capability  for  testing 
{  G r i dZoneGene ra t i on  procedures  for  TRAILBLAZER. 

< 


TYPE 

DesreesReal 

ZoneRanSe 

Letters 

VAR 

Lonsi  tude 
Latitude 
G  r i dZoneNumbe r 
GridZoneLetter 


Real) 

1  .  .60? 

'  A '  .  .  '! 


Dee  reesReal? 
DesreesReal 5 
ZoneRanSe? 
Letters? 


PROCEDURE  TBGridZoneGeneration 


<-C0UT> 
<0UT> 
•CIN  > 
■CIN  > 


VAR 

VAR 


Lone i tude 
Latitude 
G  r  i  dZoneNumbe r 
GridZoneLetter 


Dee  reesReal > 
DesreesReal ? 
ZoneRanSe ? 
Letters) ?EXTERN? 


C  PROCEDURE  TBGridZoneGene  ration  models  t'ne  TRAILBLAZER  coriversion 
-C  of  SeoSraehic  coordinates  to  Universal  Transverse  Mercator  > 


-f.  (UTM)  coordinates  erid  zone  desienator  number r letter »  > 

r  ■> 

•«.  / 

BEGIN 

.  -c  y 

URITELN  ('  ENTER  Longitude')? 

READLN  (  Lor.a i tude  )? 

URITELN  <'  ENTER  Latitude')? 

READLN  (Latitude  )? 

TBGridZoneGeneration  <  C0UT>  Lons l tude > L at i tude r 

■CIN  >  Sr idaonenumbe r * S r i dzone 1 e t te r ) ? 

URITELN? 

URITELN  ('  G  r  i dZoneNumbe  r  is  ' > G r idZwneNumber  )> 

URITELN? 

URITELN  ('  GridZoneLetter  is  '  t G r idZoneLe t te r  )? 

< 

END ,  -C  of  PROGRAM  DRVGZTB  > 


JO  /<? 


•immArnu 


rS  rS 


100 

■C 

> 

200 

MODULE  TBGZDG  <  INPUT .OUTPUT  )} 

300 

< 

I 

400 

TYPE 

500 

DeS  reesReal 

REAL  i 

600 

ZorieRanSe 

-  1  * . 60  » 

700 

Letter*  = 

'  A  '  .  .  '  I  '  i 

300 

< 

n 

r 

900 

PROCEDURE  TBGridZoneGe no  ration 

1000 

(<IN  3-  Lunditude 

*  Decrees Real. 

1100 

■CIN  y  La  t  i  uuof 

•  Did reesRee ) » 

1200 

•C0IJT3-  VAR  G  r  i  dZoneNuuiue  r  J  Zinn.‘R<jn\ite  i 

1300 

■C0UT>  VAR  G  1-iuZi.iliel.ei.Lb-  r  5  Le  t  It;  rs-  )  i 

1400 

•c 

■» 

y 

1500 

-C 

PROCEDURE  TBG  r  iuZoiieGene  r  a  L  i  on 

Models  tin.'  TRiiILBLAZER  eunver 

\>  i  u  1 1  y 

1600 

■C 

o  T  dBud  r 'it'ii i i:  coordinate*  Lu  Ui 

livar-jitl  T  r  an'.* verse  Hwreetur 

j- 

1700 

< 

(UTM)  Cuul-diiidtei  -  arid  cunt.'  de  s  i  iiiia  lu  r  uuuii.it.' r  r  letter* 

s 

1300 

■C 

y 

1900 

< 

Documental  i  un  ubeu  Hit*  GP2UM  ^ui'^  rus  i*c!m  u  L u  •  20  Tifl*  Si 

y 

2000 

■C 

"row  LIib  listings  provided  Pur 

the  TRA  11 . BLAZER  system. 

r 

2100 

-C 

s 

2200 

< 

PROCEDURE  TBGridZune Gene  rat  i on 

is  rePe renced  ha* 

2300 

< 

PROGRAM  DRVGZTB 

> 

2400 

< 

y 

2500 

■c 

PROCEDURE  TBGri JZoiit'Gt'iMfi'ciliuii 

mokes  no  re  Pe  re  no es  ♦ 

y 

2600 

< 

y 

2700 

< 

This  ppoceduit;  DOES  NOT  ft*  iTuim  any  utlc:  vul  i  ui  uv  cheeks 

y 

2800 

•c 

that  if  e  mil.  x  1  ici  tin  seec  i  P  ied  in  the  ul.iuri  thin 

y 

2900 

■c 

description*  Tiiis  is  to  allow  Hie  oi  liuri  uiw  Pee  to  re1..  In  1  *  f 

y 

3000 

■c 

rot' resen  ted  inure  el  our  Is* 

•s 

y 

3100 

•c 

> 

3200 

TYPE 

3300 

Letters  -  '  A  '  .  .  '  Z  '  » 

3400 

Inde>;Rensfe 

1..24J 

3500 

< 

■> 

y 

3600 

VAR 

3700 

Gr  idZoneL  trLis  t 

ARRiViT  1 .  .243  OF  LETTERS  » 

3300 

GridZonelndex 

IndewR.lnbe  $ 

3900 

■c 

y 

4000 

BEGIN 

4100 

< 

Initial  ice  allowable  ch* ree to i  <>  nmm 

y 

4200 

■c 

y 

4300 

G r idZoneL trL is l  TIT 

*  '  A  '  1 

4400 

G r i dZuneLl rL is l  C23 

-  '  B  '  i 

4500 

Gr  idZoneL  tr-List  C33 

- '  C  '  » 

4600 

GridZoneLlrL isl  C  4  3 

-  '  D  '  i 

4700 

Gr idZoneL LrL is t  T53 

=  ■£  ’  i 

4800 

G  r  idZoneL  t  i  l.  L * L  C6] 

-  'F '  j 

4900 

Gr idZorieLt i-List  C73 

-  '  G  '  ? 

5000 

Gr  idZoneL LrList  Co] 

- '  H  '  » 

5100 

GridZoneLtrLi  si  C93 

=  '  J  '  » 

5200 

G r IdZoneL L rL i > L  CIO] 

-  '  K  '  » 

5300 

GridZuneLLrList  C 3  1  ] 

-  '  L  '  i 

5400 

Gr idZoneL LrL 1st  C12] 

-- 

5500 

G  r  idZoneL  L  i-L  is  L  C13] 

-  '  N  '  » 

5600 

G r idZoneL t rLi st  Cl 4] 

-  •?’  i 

5700 

Gr idZoneL LrList  C153 

-'0' 1 

5800 

Gr  idZoneL  l rl. i s  t  C163 

-  '  R  '  i 

5900 

G r idZoneL L rL i s L  C173 

■=  '  S  '  J 

6000 

Gr  idZoneL  L  rl.  is  t  CIS] 

-  '  T  '  » 

6100 

G r  i uZorieL t rL i  s t  C 1 9 ] 

-  '  U  '  J 

6200 

Gr idZoneLli List  C20] 

-'V'  1 

6300 

GridZoneLti L ivt  C21] 

=  ’  W  '  » 

W 

O 

W 

O 

rS  rS 


6400  G r i iJZuueL t iL i-> l  C223  S-'X'y 

6500  G  r  i  ijluneL  t  pL  i  t,  L  C233  J-'V'f 

6600  G  i" i ijZoiiet- i, r L  i. •.>  1,  C24  3  J-'Z'» 

6700  ■£ 

6300  G i'  i i_t Zorn? Mi p  •-  TRIJNC  ( 31  ♦  0+  ( Lun-i  L  L  u  •.!*.■  /  6  .  0  )  )  5 

6900  G  r  i  dZone 1 1 idew  <  -•  TRUNC  ( 13  •  0+  (  L  <.=  i.  i  luuc/3  .  0  )  )  i 

7000  G  i"  i'iZuneLij  L  Lm-  I-  ur  LiiZuuuL  L  rl.  i.  •jLCGciuZuu<jTmit>,--;3  i 

7100  -C 

7200  END*  -C  >  j  T  PROCEDURE  TBGri'17oii*.jG>iiiL!  .*«.  ■.ion 
7300  < 

7400  END.  -C  >j?  .MODULE  GRGZTB 
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PROGRAM  DRVGZHG  (  INPUT » OUTPUT  >» 

< 

■C  F'ROGRAN  DRVGZHG  y  rov  i  dev.  <j  test  uriver 
•C  G  r  idZoiiuGeiic;  i'u  timi  i-  ruLii'iu  Pur  MAGIIC. 
< 

TYPE 


I  IJ  t  I- 1  .*  '>  l-  l  I  IvJ 


OAR 


R  3  J  i  alls 
7.  ij  1 1  e  R  j  1 1 2  e 
Letters 

V 

Loriai  t  ij  <J  e 
L  a  L  i  L  ij  d  e 
G  r  i  u Z  u i  ■  eN o m L> e  r 
G  r  i  ij’Zui  ieLe  tier 


F;  r.  A  L  i 
1  .  . 60  > 

*  A  “  .  •  '  Z  '  f 

R  a  u  i  a  1 1  s  * 

Radi  onS  r 
Z  u  1 1  tj  R  a  I  r_<  e  f 
1.  Ljllt;raf 


PROCEDURE  HGG  r  L  dZui  leGeiis  r  a  t  i  on 
(  0  0  U  T  !r  L  o  1 1  si  a  t  o  u  e 

■C  0  U  T  y  Latitude 

-C I N  y  VAR  G  i  j  dZ uiitRumLn-  r 
■CIM  y  VAR  G  r  l  dZui  ifL.;  1 1& i 


Rj'jjaiiai 

Rjdi  jiiS» 

Z  ui  i  t  R  ai  l  a  e  f 
Letters  )  5  E  X  T  E  R  N  > 


T  PROCEDURE  MoG  r  i  dZuiruGenu  ra  t  i  uu  models  tiiw  MAGIIC  '.uiiva  is  iuii 
■C  of  Sleo3reeliic  uuu  rdi  i  laLts  tu  Unjvcrvx.il  T  rails  v«  rse  riu  rviiiui' 

■C  (UTM)  vmi  I'l.iiiiatc  j  drid  jhiiiu  dua  i.  dn.j  tu  r  iiuiMi.it.' r  u  ’.utter. 

r 

A. 

BEGIN 

0 

WRITELM  ('  ENTER  Luir.ii  tudo  '  )  i 
READLN  (  l.ui is i  tude  )> 

WRITELM  ('  ENTER  Latitude')? 

READLN  (Latitude  ) ? 

MGG  r  i i.iZuiieGene  r  a  L  i  on  (  \0UT>  Loud  L  tude  i  La  t  i  tude  t 

<2N  0  f>  r  i  dZul  it  NuiuLie  r  r  G  r  i  d7ui  iifl.u  t  te 

WRITELM? 

URITELN  ('  G  r  .i  dZui  i eNuiiiue  r  is  '  r  C  r  i  dZoni.  Numi.it.' r  )» 

writeln; 

WRITELM  (  '  GridZuueLetter  j 


iurji.u.ulit'I.L'u,t'r  )  1 


END.  (  uf  PROGRAM  DRVGZMG  > 


?  C  -2  Z 
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300 
400 
500 
300 
700 
300 
700 
1000 
1100 
1200 
1300 
1400 
1500 
1600 
17  00 
1300 
1900 
2000 
2  1  0  0 
2200 
2300 
2400 
2500 
2600 
27  00 
2300 
2700 
3000 
3100 
3200 
3300 
3400 
3500 
3600 
3700 
3800 
39  00 
4000 
4100 
4200 
4300 
4400 
4500 
4600 
4700 
4800 
4900 
5000 
5100 
5200 
5300 
5400 
5500 
5600 
5700 
3800 
5900 
6000 
6100 


< 


MODULE  MGGZBG  (  INPUT* OUTPUT  )i 
0 

TYPE 

Radians 
2  o  n  e  R  a  n  3  e 
Letters 


real; 

i .  .60; 

'  A  '  «  .  '  Z 


PROCEDURE  MGGridZorieGe deration 
(<TN  >  Longitude 

TIN  >  Latitude 

-C0UT>  VAR  G  r  i  dZoneNumbe  r 


T 0 U T >  VAR  GndZoneLetter 


R  a  d  i  a  n  s  ? 

R  a  d  i  a  n  s  » 
ZoneRange? 
Letters)  ? 


•C  PROCEDURE  ?'GG  r  i  dZoneGene  ration  models  the  MA6IIC  conversion 
\  of  geographic  coordinates  to  Universal  Transverse  Mercator 
-C  (UTM)  coordinates  -  arid  cone  designator  number  £  letter. 
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•c 
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< 


Documentation  used  was  source  code  listings  from  tlie  MA6IIC  s 
d  o  c  u  in  e  n  t  CGI  03100  A  d  t  d  ♦  23  Get  l?78fpar*3*2*?0f  p5  ♦  1  67  «  }• 

y 

PROCEDURE  MGG  r  l  dZoneGene  rat  i  on  is  rePerenCw'd  by  l  > 

PROGRAM  DRVGZMG  > 

■> 

y 

PROCEDURE  MGGridZoneGeneration  m a k e s  no  references.  > 

y 

This  procedure  DOES  NOT  perform  ana  data  validity  checks  > 

that  are  not  explicitly  specified  iri  the  algorithm  > 

description,  This  is  to  allow  the  algorithm  features  to  he  > 
represented  more  clearly. 

CONST 

Pi  =  3.1415926? 

TYPE 

IndexRange  =  1..26? 

VAR 

G  r  i  d  Z  o  n  e  I  n  d  e  x  ;  I  n  d  e  a  R  a  n  g  e  ? 

BEGIN 

Calculate  the  grid  cone  number 

STATED  LONGITUDE  RANGE  IS  -lSG-PLONGITUDEPiSO  IN  DEGREES 

Gri dZoneNumbe r  :=  TRUNC<  <  { 1 30 . 0/P  i  )  *Lor.g  i  tude  + 1 30 . 0  )  /6 . 0  ) 
NO  compensation  for  wrap  around  of  grid  zone  numbers 
Deter  mi  lie  the  grid  cone  letter 

Since  no  details  are  provided  in  the  referenced  documentation 
it  is  ASSUMED  that  we  know  how  to  assign  A  or  B  fur  latitudes 
eoual  to  or  over  34  degrees  North  and  Y  or  Z  for  latitudes 
eoual  to  or  over  SO  degrees  South.  > 

y 

TRUNCATION  to  integer  is  ASSUMED  since  it  is  necessary  at  tiiiy> 
point  in  order  to  use  the  Gr idZonelnde:;  as  a  pointer.  > 

> 

STATED  LATITUDE  RANGE  IS  30-PLATITUDE-:  -84  IN  DEGREES  > 

Jo  ’  n  j  , 


ft  ✓S  r't  r't  r S  rS  r't  rS  r*t  rS  -+•  r't  rS  rS  ft  rS  ft  rt 


OiUU 

X 

6300 

Gr  idZouflnde;;  t=  TRUNC  ((<lS0.0/F'i>*Latitude  +  30.0) 

CO 

o 

6400 

r 

“S. 

y 

6  30  0 

\ 

Commute  midrange  arid  zone  letters 

j 

6600 

< 

y 

6700 

IF  G n dZone Index  <=  5 

6800 

THEN  G  r  i  dZoneLe  t  te  r  1=  CHR  (  G  r  i dZone I nde;; 

+  ORD! 

'  C  '  )  )  * 

6  900 

\ 

Here  we  Handle  the  'i‘  which  is  riut  used 

y 

70  00 

< 

y 

7100 

IF  (  G r i  d Z o n e I n d e x  >  =  6  >  AND 

7200 

(  G  r  i  d  Z  o  ri  e  I  n  d  e ;  ;  <  =  10  ) 

7  300 

THEN  GridZoneLetter  t  =  CHR(GridZoneIndex 

+  QRD  ( 

'  C  '  )  + 1  )  5 

7  400 

< 

X 

7300 

c 

Here  we  handle  the  'O'  which  is  not  used 

X 

7  60  0 

■C 

> 

7700 

IF  (  GridZonelndex  >=  11  )  AND 

7800 

(  G r i dZone Index  <=  19  ) 

7900 

THEN  GridZoneLetter  5=  CHR ( G r idZoue Index 

+  OR  D  ( 

'C' )+2) * 

3000 

r 

\ 

s 

3100 

■c 

The  rest  of  the  G r i dZone I ndex  are  biased  off  by 

ORD  (  ' 

C  '  >  > 

820  0 

■c 

y 

3300 

IF  Gr idZone Index  >  19 

8400 

THEN  GridZoneLetter  1=  CHR < G r i dZone I ndex > 

• 

3500 

-c 

y 

3600 

r 

\ 

Assign  Y  or  Z  to  the  North  Folar  Zone  according 

as  Western  > 

3700 

-c 

or  Easte  r  n  Hemisphere*  respectively. 

y 

3800 

•c 

y 

3900 

IF  Latitude#(130.0/Pi)  > =  34.0 

9000 

THEN  IF  Longitude*! 180.0/Pi )  <  0.0 

9100 

THEN  GridZoneLetter  1=  ' Y ' 

9200 

ELSE  GridZoneLetter  t  ■=  'Z'* 

9300 

-c 

y 

9400 

< 

Assign  A  or  B  to  the  South  Polar  Zone  according 

as  Western  > 

9500 

•c 

or  Eastern  Hemisphere*  respectively. 

y 

9600 

-c 

y 

9700 

IF  .atitude*! 180.0/Pi >  <=  -30.0 

9300 

THEN  IF  Longitude*! 130.0/Pi)  <  0.0 

9900 

THEN  GridZoneLetter  1=  'A' 

10000 

ELSE  GridZoneLetter  :=  'B'» 

10100 

< 

■%. 

j 

10200 

NO  correction  for  the  four  irregular  zones  - 

y 

10300 

-c 

32X  *  34X  *  and  36X  do  riot  exist 

y 

10400 

-c 

31V  is  truncated 

> 

10500 

■c 

\ 

10600 

END*  -C  of  PROCEDURE  M  G  G  r  i  D  Z  o  n  e  G  e  n  e  r  a  t  i  o  n 

*\. 

j 

10700 

S 

y 

10800 

END.  -C  of  MODULE  MGGZDG 

■N 

s 

J(5  y 


rS  fS  rS  rS  rS  rS  rS  <S  rS  rS  rS  fS  rS  rS  fS  r\  rS  /N  r\  rS  r'"»  rS  rS  rS 


100 
200 
.300 
400 
500 
600 
700 
800 
900 
1000 
1100 
1200 
1300 
14  00 
1500 
1600 
1700 
1300 
1900 
2000 
2100 
2200 
2300 
24  00 
2500 
2600 
2700 
2800 
2900 
3000 
3100 
3200 
3300 
3400 
3500 
3600 
370  0 
3800 
3900 
4  000 
4100 
4200 
4300 
4400 
4500 
4600 


PROGRAM  DRVGZBT  (  I NF'UT » OUTPUT  )? 

< 

-Z  PROGRAM  DRVGZBT  provides  a  test  driver  cape bilibs  for  testing 
\  GridZoneGeneralion  procedures  for  BETA. 

-C 

TYPE 

Radians  —REAL? 

ZoneRanSe  =  1..6C? 

Letters  =  ' A  '  •  •  '  Z ' ? 

VAR 

Longitude  1  Radians? 

Latitude  !  Radisns! 

G  r  i  d  Z  o  n  e  N  u  m  b  e  r  1  ZoneRanSe? 

GridZonw Letter  1  Letters? 

CenterMeridian  1  Radians? 


PROCEDURE  BTGridZoneGeneration 
( -C0UT>  Lonsitude 

-CGUT}  Latitude 

-CIN  >  VAR  G r i dZoneNumbe r 
-CIN  >  VAR  GridZoneLetter 
-CIN  >  VAR  CenterMeridian 


J  R  a  d  i  a  n  s  ? 

1  Radians? 

1  ZoneRanSe  ? 

1  Letters? 
i  Radians) 5EXTERN? 


•C  PROCEDURE  BTGridZoneGeneration  models  the  BETA  conversion 
-C  of  Seoaraphic  coordinates  to  Universal  Transverse  Mercator 
-C  (UTM)  coordinates  arid  zone  designator  number ' 1 ette r  and  the 
-C  central  meridian  of  the  rectangle. 

-C 

BEGIN 

r 

~K. 

URITELN  ('  ENTER  Lonsitude'.)? 

READLN  <  Lonsitude  )? 

URITELN  ('  ENTER  Latitude')? 

READLN  (Latitude  )? 

BTGridZoneGeneration  \  -CQUT>  Lons  i  tude  »  Lat  i  tude  » 

■CIN  >  G r i dZoneNumbe r » G r i dZoneLe t Le r 
CenterMeridian) ? 

URITELN? 

URITELN  ('  GridZoneNumber  is  '» G ri dZoneNumbe r  )? 

URITELN? 

URITELN  ('  GridZoneLetter  is  '»GridZoneLetter  )? 

URITELN? 

URITELN  < '  CenterMeridian  is  ''CenterMeridian  ) 

C 

END.  -C  of  PROGRAM  DRVGZBT  > 


3  o  - 1  y 


fS  rS  rS  M 


100 

200 

300 

400 

500 

600 

700 

800 

900 

1000 

1100 

1200 

1300 

1400 

1500 

1600 

1700 

1300 

1900 

2000 

2100 

2200 

2300 

2400 

2500 

2600 

2700 

2800 

2900 

3000 

3100 

3200 

3300 

3400 

3500 

3600 

3700 

3800 

3900 

4000 

4100 

4200 

4300 

4400 

4500 

4600 

4700 

4800 

4900 

5000 

5100 

5200 

5300 

5400 

5500 

5600 

5700 

5800 

5900 

6000 

6100 

6200 

6300 


MODULE  BTGZDG  (  INPUT. OUTPUT  )? 

< 

TYPE 

Radian* 

Zuiiv/Runae 

Letter* 

PROCEDURE  BTGr iuZuneGene  ra  Liun 


REAL. 

1 .  •  6  0  » 

A  '  »  .  r. 


(-CIN  y  Luuiii  Lude 

■C I N  >  Latitude 

-COUT>  VAR  G  r  L  dZwi leNumue  r 
<OUT>  VAR  G  r  idZutieLe  L  Le  r 
<OUT>  VAR  Cun to rMe r i d i an 


R  a  d  L  a  1 1  s  ? 
Radian*. 
Zone Rand 6  ? 

Lo  !■  i  t.  rs.  ? 
Rad  i  si i*  )  » 


-C  PROCEDURE  8TGr  LdZunoGenerstion  mode  1  s  the  BETA  c  ui  i  Vi-;  r  *  i  Ofi 
•C  of  deodraphic  cun  r  d  i  na  Le  *  tu  Universal  Transverse  Me  rea  It*  r 
■C  (UTM)  CuO  ruins  Lws  -  if  r  i  d  tunu  das  i  di  is  lu  r  iiumiLici  .  lottui  .  and 
■C  the  central  meridian* 


■C  Documentation  used  was  source  code  listings  from  the  BETA 
■C  document  8322-43  dtd.  16  Oct  1931.  on;-; » 4 .  pii  .  2  —  474  fur  the 
■C  ADSCNU  sulip  rod  ran  and  ri.i.2  — 430  for  the  ADSCCii  suin'  rut  ram 
< 

•C  PROCEDURE  BTG  r  idZuntGene  r  a  t  i  on  is  r  c  fc  reuceu  uui 

<  PROGRAM  ORVGZBT 
■C 

{  PROCEDURE  BTGr  idZone  Gem?  rat  inn  imika*  no  re  Pe  ranees  » 

■C 

■C  This  procedure  DOES  NOT  perform  .nra  data  valid  its  cheeks 

■C  that  are  not  explicitly  spec  i  f  i  ed  in  the  aliioritiiut 

-C  deso  r  Lr  t  ion.  This  is  tu  allow  the  algorithm  features  to  be 

<  represented  mure  cl earl  a* 


<  Since  the  included  '  ZDBPR0 « COM '  is  nut  available  t 
-C  this  time,  we  assume  implicit  tueinii  in  Lite  source 
■C  code. 


to  Ub  a  u 

a  FORTRAN 


CONST 

Pi 

TYPE 

Letters 
X  ndexRansle 


Gr  i dZoneL Lr List 
G  r i dZunel nd ex 


=  3.1413926? 


A  '  *  .  c 
1 .  .  2  4  ? 


ARRAY: 1.. 24 3  OF  LETTERS? 
1 1 idexRsi  ixie  ? 


BEGIN 

Initialize  allowable  characters  orrai 


GridZorieLtrL  it.  t  Ml 
G ri  dZoneL. t rl.i  * t  C21 
G r i dZoneL t rL i s t  131 
0  r  i  dZoneL  L  rl.  i  s  L  C41 
GridZorieLtrL  i  at  C31 
G  r  i  dZoneL  t  rl.  i  s  t  Col 
Gr iuZorieL t rList  C71 
G  r  i  dZoneL  t  rl.  i  s  t  C81 
Gr idZoneL t rLis t  C?1 
Gr  idZoneLlrl.isL  '10! 


GridZorieLtrL  is t  Clll 


JO  -X(c 


S'  S'  s 


6400 

G  r  idZunel.  t  r L  i a  t 

C  123 

- '  M  '  5 

6500 

GridZoneLLrt  i  a  t 

C  1  33 

-  '  N  '  5 

6600 

G  r  id Zui ie(.  L  r L  L  a  t 

C  1 4  3 

-'P'5 

6700 

GridZuneLtcLial 

C  133 

-  '  Q  '  5 

6800 

G  r  idZuneL  t  rL  i  s  t 

C  1 6  3 

-  'R'  5 

6900 

G r i dZoneL  t  rL 1s  t 

C  3  7  3 

=  '  3 ' 5 

7000 

G  r  i dZune  L  t  r  L  1st 

E  1 3  3 

• '  T  '  5 

7100 

GridZoneLtrL  ist 

Cl  93 

-  '  U  '  * 

7200 

G  r  i  d/oneLt rL ist 

C20  3 

'V'  5 

7300 

G  r  i  dZone L  t  rL  i  a 1 

C23  3 

^'W'i 

7400 

G  r  i  dZoneL  t  r  1.  1st 

C  22  3 

-'X'  5 

7500 

GridZuneLLrList 

C  23  3 

—  '  Y  '  5 

7600 

G  r  i  dZoneL  t  r L  j.  s  t 

C  24  3 

-'Z'  5 

7700 

-c 

> 

7800 

< 

Calculate  the  'Jrid  z 

one  number 

J* 

7900 

< 

> 

3000 

G  r  i  dZunoNumbu  r 

i  -  <  TRUMC  <  (  ( 1 3000 . 0/F’i  )  £Lon-.<  i  lude ) 

8100 

+18600.0))  DIV  6005 

8200 

< 

J 

8300 

< 

Compensate  Fur  wre*- 

a  round  uf  arid  zone  numbers  > 

3400 

•c 

8500 

IF  G  r  i dZuneMumi.it/ 

r  >  60 

3600 

THEN  GridZonc 

Numb*;  r  •  =  GridZuneNumber_60 * 

8700 

8300 

8900 

9000 

9100 

9200 

9300 

9400 

9500 

9600 

9700 

9800 

9900 

10000 

10100 

10200 

10300 

10400 

10500 

10600 

10700 

10800 

10900 

11000 

11100 

11200 

11300 

11400 

11300 

11600 

11700 

11800 

11900 

12000 

12100 

12200 

12300 

12400 

12500 

12600 


■C 

< 

< 

< 

< 

-c 

■c 

■c 

-c 

r 

■C 


< 

< 

C 


IF  G  i-  i  dZuneNumbe  r  <  1 

THEM  G r  l  uZuneNumba r  «-  G  r idZoueNumbe  r  +  60  5 

Do tw i‘ im  1  m.1  tin;  dr  id  zomi  letter 

G  i"  IdZoiio  IiiiJe;-;  i  -  TRUMC  (  ( I.  >i  t  j.  tuu*;$  (  180*0/ Pi  )+104.0)/8. 

f 

Test  f  u  i‘  ■Hid  lock  uut  Nurtn  color  Zones 
IF  G  »•  i uZune Index  >  22 

THEN  G  ridZouelude*  J  22* 

NOTE  that  no  such  lest  is,  needed  Fur  L-  i  i  c-  South  Pul  or  Zulu.' 
because  the  olUoriim  limit  woo  :iivon  -su  SO  Sooth 


G  r idZonol.o  t  It;  r 


l*  r  i  d  Zonal.  t  i‘L  i  s  LF  G  r  j.  uZoi  le  I  not;;;  I  * 


Cor  rue t  For  the  Four  i  r  cedilla r  aunt's  - 
32X*34X*  and  36X  do  nut  exist 
31 V  is  truncated 

T rune ate  dr  id  <:ont)  31V 

IF  < G r i dZune l ndex  -  20)  AND 
< < G r i dZuneNumbe r  -31)  AND 
(l.ondituue  >-  3  ♦  0#  ( P  i/130  *  0 ) ) ) 

THEN  G  r  i  dZuneNumbe  »•  5  32  5 

Correct  for  arid  /.ones  32X*34X*  end  3oX 

IF  <  G  r  i  dZune  Index  -22)  AND 
( G r i dZuneNumbe r  -32) 

THEN  IF  Lunsitude  >-  9  *  0* < P i / 1 80 . 0 ) 
THEM  G  r  i  dZoneNumb*;  r  »-  33 
ELSE  G  r  j  dZuneNumbe  r  J  31* 


IF  <  G  r  idZune Index  -  22)  AND 
< G r i dZuneNumbe r  -  34) 

Jo-,’7 


V»  s>  S'  S'  V  S'  V  S'  V  S'  S'  S'  S'  S'  S'  S'  S' 


12700 

12300 

12900 

13000 

13100 

< 

13200 

IF 

13300 

13400 

13500 

13600 

13700 

< 

13300 

Cei  i 

13900 

•C 

14000 

emi  > 

-C  o  f 

14100 

■C 

14200 

END . 

■C  of 

THEN  IF  Lunaituu'v  >-  21 . 0*  < F i / 1  SO  .  0  ) 

THEN  6  r  i uZuneNuiiibe  r  t-  35 
ELSE  C5  r  i liZuiioNuiuiji.'  r  • 33 » 


(GriuZune Number  -  36) 

THEN  IF  LoiWLluuu  >-  33 . 0$  ( P  i  / 1 30 . 0  > 

THEN  G  r  i  uZuneNuuiht' r  J-  37 
ELSE  G  r  i  iJZoi  ieNumbe  r  t  -  35? 

Cent,  e i' Me r  .i. d i.  tm  J  •-  <  6*G  r  i  liZorinMumbe  r  - 1 83 )  #  ( P  i.  / 130  »  0  ) 


> 

> 

> 


3  o  -if 


100 

200 

300 

400 

300 

700 

800 

900 

1000 

1100 

1200 

1300 

1400 

1500 

1600 

1700 

1750 

1300 

1900 

2000 

2100 

2200 

2300 

2400 

2500 

2603 

2604 
2606 
2608 
2900 
3200 
3300 
3400 
3500 
3600 
3700 
3300 
3900 
3950 
4000 
4100 
4200 
4300 
4400 
4500 
4550 
4600 


•  Du  J  rytf'jR'.'  a  1  ? 

!  flwsi  reel*  Real  5 
'  Zui  m  R a  mi e  ? 

*  Lt-LLd-v.  f 

t  Dual' «  a  1  ii  L  a  :iu  r  3  5  E  a  T  E  R  N  ? 


PROGRAM  DRVGZGR  (  INPUT » OUTPUT  3 ? 

■C  > 

•C  PROGRAM  DRVGZGR  pruviues-  a  LeuL  ui-ivt'i'  curei*  1  1  i  La  Pur  LeuLimi  > 
■C  G  r  i  dZone  G<#i  ia  r  a  L  i  ui  i  Pur  GUARDRAIL.  s 

<  > 

TYPE 

DedreeaReal  =  Rea-3  ? 

Ded I'eiiiilri Le de r  -  INTEGER? 

ZorieRande  =  1..60? 

L  a  tiers  -  '  A ' . • ' Z ' ? 

VAR 

Longitude  J  Iit.-urefaRi  cd  ! 

La l i Lude  t  Dad reeaRea 1 ? 

G  r  idZoiieNumlie  r  J  ZuneRemie ? 

G r i dZuneLe L Le r  t  Looters? 

ceri  Le  rme  r  i  d  i  an  J  dud  reet.  i  ii  Lt-de  r  ? 

<  > 

PROCEDURE  GR  G  i*  i  dZ  une  G  e  ne  r  a  L  i  un 

( {OUT}-  Lund i Lude  «  DudruesReul  ? 

-C0UT>  La Li Lude  *  Dud reusRue 3  ? 

■CIN  >  VAR  Gr  LdZoiieNumuer  '  ZuiiuRumie ? 

■CIN  }■  VAR  G  r idZoneL  e  L  La  i  *  Ll  L  Le  rv»  ? 

■CIN  s  VAR  Ciin  Le  r  Me  r  i  d  i  an  ?  Dud  ieav>I  ii  Lada  r  3?  EXTERN  ? 

<  > 

■C  PROCEDURE  GRG  r  L dZuimGene  im L  l uu  Mudclu  Lite  GUARDRAIL  i.ui ivu  r  ->  i un  j 
■C  of  aeod  i-ap-hic  CuUrdineLet,  Lu  Uii  i  vt  i-ua  1  Traii4.vtM-t.fe'  Me  n-t  Lo  <•  > 

C  <UTM)  i.-oo  rd  i  i  id  Leu  -if  id  <'.on>.>  idne  Lu  r  mi#o«  r  f  la  L  Le  r  .mu  Lite 

■C  central  meridian  uP  Lite  rttLciialv.  > 

C  > 

BEGIN 

-C  > 

URITELN  ('  ENTER  Luudi tude ' 3 ? 

READLN  (  Luna j Lude  3  5 
WRITEI.M  ('  ENTER  Latitude  '  3  ? 

READLN  (Latitude  3? 

GRG  i* i dZoueOene  ra  L i  on  (  -C0UT3-  Lund  i  Lude  f  La  L  i  Lude  » 

■CIN  y  d  v  i  d.c:uiienuiiilie  r  r  r  i  d.:une  1  v  L  Lu  r  f 
Can  La  n»e  r  t  d  i.  an  3  ? 

URITELN? 

URITELN  ('  G r  i dZuueNuMbe r  i».»  '  t  G  r  i  dZulieNuuiUe  r  3? 

URITELN? 

URITELN  ('  G l  idZuneLe  L  Le  r  i  i>  '  »  G  r  idZunel.  e  L  Li:  r  3? 

URITELN? 

URITELN  ('  CenLu  r Me  r  i  d  i  all  j  <>  '  F  Cel i  Le  i-Mu  r  i  u i  ail  3 
<  3 

END.  -C  of  PROGRAM  DRVGZGR  > 


3o 'X  ^ 


100 

200 

300 

400 

500 

600 

700 

300 

900 

1000 

1100 

1200 

1300 

1400 

1500 

1600 

1700 

1800 

1900 

2000 

2100 

2200 

2300 

2400 

2500 

2600 

2700 

2300 

2900 

3000 

3100 

3200 

3300 

3400 

3500 

3600 

3700 

3300 

3900 

4000 

4100 

4200 


MODULE  GRGZDG  <  INPUT, OUTPUT  >1 

< 

TYPE 

DedreesReal 
Ded  ree-aln  Cede  r 
ZoneRande 
La tiers 

< 

PROCEDURE  GRG  r  i  dZuneGen:-1  r  a  l  i  on 


REAL! 

integer; 

i . .60; 

'  A  '  .  .  '  Z  ' 


(  \  I N  >  Lundi lude 

•CIN  >  Latitude 

•£0UT>  VAR  G  riuZuiieNuiiii.it- 1 
■CGUT3-  VAR  GridZuneL filer 
<0UT>  VAR  Ccn le rMe r i d i an 


Ded  r eesRea  1  i 

Oe:i  1'L‘ii’iRii-j  1  J 

Z  u  1 1  e  R  a  1 1  d  v  i 
L  e tiers; 

Ded  iff  S  I  ii  Cede  i- ) 


< 

< 

< 

< 

< 

< 

< 

< 

< 

< 

< 

•C 

-c 

< 


PROCEDURE  GRGr idZuneGener  lion  mu',  it' Is  the  GUARDRAIL  eunve  rsi 
of  S'-'O-J  r  aid  i  i  e  cuo  r  d  i  lie  ley  to  Univui">iil  T  r  Jiijvi;  r  \>e  Me  rue  lu  r 
(UTM)  coordinates  -  dlid  acme  desidnalur  miiiutTi  letter,  and 
tile  central  meridian. 


PROCEDURE  GRGriuZune Gene  rat  ion 
PROGRAM  DRVGZGR 


i  a  I' e  l  ti  l'll i; f  IJ  US. 


PROCEDURE  GRG  r  i  uZuneGeiie  i  a  t  c  oi i  makes  nu 


re  i  e  re  net1  a  * 


s 

Cl  I  l}- 

y 


j 

s 

s 


This  procedure  DOES  NOT  reiTum  any  data  validity  cheek1, 
that  are  not  explicitly  spec i P iod  in  tin;  alioritimi 
description.  Tiiis  is  to  allow  the  aldoritiiM  "i.'i:turt!)  to 
ree reseri tod  more  clearly. 


> 

}• 

> 

> 


TYPE 

Letters 

IndexRenSe 

VAR 

GridZoueLtrList 
G  r  i dZune  Index 

BEGIN 

Initialize  allowable 


'  A  '  .  »  '  i 
i . .  2i ; 


ARRAY! 1 . .241 

1 1 idesRai uie  f 


op  letters; 


ena rac ec rs  or 


4300 

G  r idZuneL l  i-List 

rn 

- 

'A  ' 

4400 

G  r  iuZonuL  t  i-l.is  t 

C  21 

- 

'3' 

4500 

G r idZuneL trList 

T33 

-= 

'C' 

4600 

6  r idZuneL t  rL L a  t 

C43 

-- 

'0' 

4700 

G  r  idZuneL 1  rL  i  s  t 

rs: 

- 

'E  ' 

4300 

G  r  L  i.iZornit.  t  rL  Lai. 

C  6  3 

- 

4900 

G  r idZuneL trList 

C  7  3 

- 

'G  ' 

5000 

G  r  i  d  Z  u  n  e  L.  t  r  L  Lit 

C33 

■  - 

'  H  ' 

5100 

G  r  idZonel.  t  rL  i  s 1 

C  9  3 

- 

'  J' 

5200 

GridZuneL  t rL i s t 

C  1 0  3 

■- 

'  K  ' 

5300 

GridZuneL trList 

cm 

- 

'L  ' 

3400 

G  r  i  dZui  lel.  t  rL  L  •>  l 

ri23 

- 

'  M  ' 

5500 

GridZoueLtrList 

Cl  33 

■3 

'H' 

5600 

G  r  i dZuuel.  t  rL  Let 

C  1 4  3 

- 

5700 

GridZuneL trList 

C  1  33 

'  R  ' 

5800 

Or  i  dZui  i«L.  t  rL  L  a  t 

C  1 6  3 

■  - 

'R  ' 

5900 

G  r idZuneL t rL i s  t 

Cl  73 

- 

'S' 

6000 

G  r  i  dZuneL  t  rL  i  a  l 

C  133 

-- 

'T' 

6100 

GridZuneLtrList 

C  1 9  3 

- 

'  U  ' 

6200 

0  r i dZnueL t rL  L  a  t 

C  20  3 

'  V  ' 

6300 

GridZuneLtrList 

C21]JCdo 

a 

'VI ' 

ay 


<4400 

G  i"  i  u  ZorneL  t  i-L  i  L 

c  22 1 

6500 
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